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A  discrete-event  simulation  of  a  stochastic  process,  a 
machine-repairman  model,  has  been  programmed  on  the  IBM 
Personal  Computer.  The  model  consists  of  three  helicop- 
ters, of  which  two  are  in  service  and  one  is  in  cold 
standby,  with  an  option  of  one  or  two  repairmen. 

The  program  output  is  a  graphics  display  containing  a 
system  state-versus-time  graph,  a  table  of  statistics,  and 
animated  figures  that  illustrate  the  current  state  of  the 
system.  The  program  user  can  directly  observe  the  dynamics 
of  the  model  as  the  fixed-increment,  simulation  clock 
advances.  The  user  has  the  option  of  changing  the  follow- 
ing model  parameters:  helicopter  failure  rate,  repairman 
service  rate,  and  the  number  of  repairmen  to  employ. 
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I.   INTRODUCTION 

During  recent  years,  the  versatile  microcomputer  has 
found  a  niche  in  many  of  the  professional  discliplines 
including  education.  For  difficult  and  complex  subjects, 
educators  are  exploring  the  possibilities  of  the  microcom- 
puter as  a  teaching  aid  in  addition  to  its  more  traditional 
role  as  a  computing  machine.  Some  universities  now  use 
microcomputer  graphics  to  illustrate  principles  that  cannot 
be  fully  explained  in  a  lecture.  [Ref.  1] 

A  widely  used  application  of  computers  is  to  imitate  or 
simulate  stochastic  processes.  This  technique  usually 
involves  only  the  numerical  evaluation  over  time  of  the 
corresponding  model  with  the  intent  of  gaining  some  under- 
standing of  the  process.  But  using  graphics  as  part  of  the 
output  in  a  simulation  provides  the  user  another 
perspective.  For  example,  observing  a  working  model--such 
as  a  birth-death  process--can  help  in  grasping  the  funda- 
mentals of  the  model  and  provide  some  insights  into  the 
process  not  easily  seen  on  the  classroom  chalkboard. 

This  paper  describes  the  development  of  a  discrete- 
event,  graphic  simulation  of  a  specific  example  of  the 
machine-repairman  model.  The  simulation  gives  the  observer 
the  impression  of  real-time  passage.  A  graph  on  the  visual 
display  shows  the  history  of  the  system's  transitions  from 


state  to  state,  and  animated  figures  illustrate  the 
dynamics  of  the  model.  The  simulation  also  exhibits  a 
table  of  the  theoretical,  and  the  current,  estimated  values 
of  quantities  commonly  used  in  evaluating  queueing  systems. 

Section  II  of  this  paper  is  a  discussion  of  the 
machine-repairman  model,  the  assumptions  upon  which  the 
model  is  based,  and  the  methods  used-  for  obtaining  the 
theoretical  and  estimated  values  of  the  quantities  of 
interest.  Section  III  covers  the  simulation's  implementa- 
tion on  the  IBM  Personal  Computer  (IBM  PC).  Appendix  A  is 
a  user's  guide  which  contains  summary  information  about  the 
machine-repairman  model  and  instructions  for  operating  the 
program.  The  program  flowchart  and  the  program  listings 
are  found  in  Appendixes  B,  C,  D,  and  E. 


II.   THE  MACHINE-REPAIRMAN  MODELS 

A.   DISCUSSION 

The  simulation  described  in  this  paper  actually  con- 
tains two  machine-repairman  models.  In  both  models  there 
are  three  machines,  represented-  by  helicopters.  The  models 
are  distinguishable  by  the  number  of  repairmen;  there  is 
one  repairman  in  the  first  and  two  repairmen  in  the  second. 

The  nominal  mission  is  to  keep  two  helicopters  flying 
at  all  times.  The  third  helicopter  is  placed  in  cold 
standby  as  a  backup  to  the  other  two.  If  available,  a 
repairman  is  assigned  immediately  to  a  helicopter  upon  its 
failure.  If  a  repairman  is  not  free,  then  the  helicopter 
joins  a  queue  to  await  repair. 

Several  assumptions  are  made  to  simplify  the  models.  In 
the  two-repairmen  model,  repairmen  are  assumed  to  be 
equally  competent  and  to  take  the  same  mean  time  to  repair 
a  helicopter.  Both  models  assume  the  repair  time  to  be 
continuous;  that  is,  repair  parts  are  always  available,  so 
there  is  no  delay  in  repairing  a  machine  other  than  the 
time  required  by  the  repairman  to  perform  the  work.  The 
repairmen  work  independently  and  do  not  assist  one  another. 
Similarly,  the  models  assume  the  helicopters  operate  inde- 
pendently and  have  the  same  service  life  with  the  same  mean 
time   to   failure.   By   assumption,   a   helicopter   is   not 
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subject  to  failure  while  in  cold  standby  since  it  is  not 
operating. 

A  machine-repairman  model  typically  assumes  its 
machines  operate  continuously  until  failure.  In  the  case 
of  aircraft,  this  assumption  does  not  reflect  reality  since 
in  actual  practice  aircraft  do  not  fly  continuously,'  but 
have  instead  periods  of  inactivity  between  missions  of 
limited  duration.  As  already  noted,  one  can  reasonably 
assume  that  an  inactive  helicopter  is  not  subject  to  fail- 
ure and  that  the  failure  process  occurs  only  when  the 
helicopter  is  active. 

The  actual  failure  rate  of  a  flying  helicopter  is  the 
number  of  failures  per  flight  hour.  The  product  of  the 
actual  failure  rate  times  the  ratio  of  the  number  of  flight 
hours  performed  to  the  number  of  calendar  hours  elapsed 
until  failure  defines  the  helicopter's  effective  failure 
rate.  Thus  the  effective  failure  rate  is  the  number  of 
failures  per  calendar  hour.  To  simplify  the  simulation, 
the  models  assume  that  the  helicopters  operate  continuously 
until  failure  with  the  effective  failure  rate. 

The  selection  of  helicopters  to  represent  the  machines 
in  these  models  allows  a  simple  graphic  design  in  depicting 
their  status.  An  operating  machine  is  seen  in  the  display 
as  a  green,  flying  helicopter;  a  helicopter  that  is  in  cold 
standby  is  also  colored  green  but  is  not  flying.   A  machine 
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that  is  down  and  awaiting  repair  is  colored  red.   When  the 
down  machine  is  in  repair,  it  is  seen  with  a  repairman. 

Figure  1  is  a  frequently  used,  pictorial  representation 
(see  [Ref.  2],  for  instance)  of  the  machine-repairman  model 
with  one  repairman.  The  dashed  lines  enclose  the  model's 
repair-queue  system  which  contains  a  single-server  queue. 
The  machines  outside  the  dashed  box  are  the  potential 
customers  for  the  repair  queue.  Since  there  are  only  three 
machines  in  the  model,  the  longest  possible  length  of  the 
queue  will  be  two. 
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Figure  1:  Machine-Repairman  Model  with  One  Repairman. 

This  model  has  a  second  queue,  the  cold-standby  queue, 
since  a  helicopter  may  have  to  wait  before  being  assigned 
to  a  mission.  The  standby  queue  can  contain  only  one  heli- 
copter, and  this  occurs  when  there  are  no  machines  in  the 
repair-queue  system. 
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The  cyclic  nature  of  the  model  is  clearly  evident  in 
the  figure.  A  helicopter's  life  cycle  consists  of  an  oper- 
ating period  followed  by,  after  a  possible  delay  in  queue, 
a  repair  period  which  is  in  turn  followed  by,  after  another 
possible  delay  in  cold  standby,  another  operating  period. 
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Figure  2:  Machine-Repairman  Model  with  Two  Repairmen. 

Figure  2  depicts  a  pictorial  representation  of  the 
machine-repairman  model  with  two  repairmen  [Ref.  3].  In 
this  system,  helicopters  waiting  in  the  repair  queue  go  to 
the  first  available  repairman  for  service.  The  maximum 
possible  length  of  the  repair  queue  is  only  one.  As  in  the 
case  of  the  one-repairman  model,  the  standby  queue  can 
contain  at  most  one  helicopter. 

The  machine-repairman  model  is  a  special  case  of  a 
birth-death  process  in  which  the  event  of  a  machine  failing 
at  time  t  constitutes  a  birth.   Let  T  be  the  operating  life 
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of  a  helicopter  until  failure  and  assume  that  T  is  exponen- 
tially distributed  with  a  failure  rate  a.  Then  the  mean 
time  to  failure  (MTTF)  of  a  helicopter  is  E(T)  =  1/a. 
Similarly,  the  completion  of  repair  on  a  machine  is  a 
death,  so  define  S  to  be  the  service  time  of  a  repairman. 
The  service  time  is  also  assumed  to  have  an  exponential 
distribution  with  a  service  rate  y.  The  mean  time  for 
repair  (MTFR)  is  then  E(S)  =  1/y. 

For  these  machine-repairman  models,  the  state  of  the 
system,  N(t)  =  n,  is  defined  as  the  number  of  helicopters 
down  at  time  t.  Since  there  are  only  three  helicopters  in 
the  models,  then  n  =  0,1,2,3. 

B.   MODEL  ESTIMATORS 

In  the  study  of  queueing  models,  there  are  four  quanti- 
ties that  are  commonly  used  to  evaluate  the  queueing  sys- 
tem's performance.  As  part  of  the  graphical  output,  the 
simulation  will  provide  the  user  a  table  of  theoretical  and 
estimated  values  of  these  four  quantities: 


L   =  the   average   number   of   machines   in   the   repair 
system, 

L   =  the   average   number   of   machines   in   the   repair 
q   queue, 

W   =  the  average  amount  of  time  a  machine  is  down, 

W   =  the  average  amount  of  time  a  machine  waits  for 
repair. 
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As  shown  in  Figures  1  and  2,  the  variable  q  is  the  time 
that  a  down  machine  waits  in  the  repair  queue.  Then 
W   =  E(q).   It  follows  that  the  average  total  time  that  a 

q 

helicopter  is  in  the  repair  system  is  the  sum  of  the 
average  waiting  time  in  queue  and  the  expected  service 
time:  W  =  W   +  E(S  )  . 

q 

Little's  formulas  show  the  relationship  of  the  four 
quantities  of  interest 

L  =  AW 
L   =  XW 

q    q 

where  X  is  defined  as  the  average  arrival  rate  of  machines 
entering  the  repair-queue  system. 

There  are  three  additional  quantities  that  will  be 
computed  and  presented  by  the  simulation.  They  are  defined 
as 

P,  =  P{down  machine  must  wait  for  repair}, 

P   =  P{up  machine  must  go  to  cold  standby}, 

Av  =  availability   =  the   proportion   of   time   that   at 
least  one  helicopter  is  not  down. 

To  calculate  all  of  eight  of  these  values,  it  is  first 
necessary  to  compute  the  system's  steady-state  probabili- 
ties, p  =  P{N  =  n},  n  =  0,1,2,3.  In  other  words,  p  is 
the  long-run  probability  that  the  system  is  in  state  n.  The 
details   of   calculating   p  ,   a   function   of   the   machine 

3     C  n 
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failure  rate  and  the  repairman  service  rate,  will  not  be 
discussed  in  this  paper.  The  interested  reader  should 
refer  to  the  bibliography  for  a  list  of  textbooks  on  queue- 
ing  theory  and  birth-death  processes.  In  addition  to  the 
limiting  probabilities,  {p  ,  n  =  0,1,2,3},  the  program 
computes  the  fraction  of  the  time  the  system  is  in  state  n, 
denoted  p  ,  as  equal  to  the  total  time  in  state  n  divided 
by  the  total  elapsed  time.  The  quantity  p  is  an  estimate 
of  pn. 

The  value  of  L,  the  theoretical  average  number  of 
machines  in  the  repair  system,  is  found  by  averaging  the 
number  of  down  machines  over  all  states: 


L  =   )  np   =  (l)p,  +  (2)p,  +  (3)p 
n4o   n       x       *       J 


The  estimate  of  L,  denoted  L,   is  computed  with  the  same 

expression  with  the  exception  that  the  estimate  p   is  sub- 
stituted for  p  . 
^n 

The  average  number  of  machines  in  the  repair  queue  L 

and  its  estimate  L  depends  in  part  upon  the  number  of 

repairmen   in   the  model.   In   the   one-repairman   case, 

machines  occupy  the  repair  queue  only  when  the  system  is  in 
state  2  or  3 ,  thus 


Lq  =  P2  +  2p3,  and 


L   =  f>~    +    2p. 
q    v2  . 
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In  the  two-repairmen  case,   a  queue  forms  only  when  the 
system  is  in  state  3,  thus 


L   =  P-,  ,  and 
q   ^3 ' 


L   =  p-  . 

q    ^3 

The  average  rate  at  which  helicopters  enter  the  repair 
queue,  \,  is  dependent  upon  the  number  of  active  machines 
subject  to  failure  in  state  n  and  the  proportion  of  time 
the  system  is  in  state  n.  There  are  two  active  helicopters 
in  states  0  and  1,  one  in  state  2,  and  none  in  state  3.  The 
combined  failure  rate  of  the  active  helicopters  in  state  n 
is  given  by 


a   =  < 


2a 

n    =    0,1 

a 

n    =    2 

0 

n    =    3. 

Then  X  is  found  by 


V 
„^0 


U    anpn   =  2ap0  +  2opl  +  a*2 


For  the  estimate  X,  the  program  divides  the  number  of  fail 
ures  that  have  occurred  in  time  t  by  t ,  that  is 

X  =  total  number  of  failures/total  elapsed  time. 


17 


The  average  waiting  time  quantities,  W  and  W   are  then 
calculated  by  using  Little's  formulas 


W  =  L/X 
Wq  =  Lg/X. 

To  compute  the  estimate  W  ,  the  program  keeps  a  record  of 
the  length  of  time  each  helicopter  spends  in  the  repair 
queue  and  uses  the  expression 


m 
W   =  (1/m)  )   q.      m  =  0,1,2, . . . 
q        i%l 


where  q.  is  the  waiting  time  in  queue  of  the  ith  helicopter 
and  m  is  the  total  number  of  helicopters  that  have  entered 
the  queue  by  time  t.   Then  W  is  easily  found  with 

W  =  W   +  E(S)  . 

q 

The  estimate  P,  of  the  probability  that  a  down  machine 
must  wait  for  repair  is  the  proportion  of  the  time  that  the 
system  spends  in  a  state  in  which  a  machine  is  down  with  no 
unemployed  repairman.  The  repairman  is  not  available  in 
the  one-repairman  model  when  the  system  is  in  state  2  or  3 : 

Pd  =  P2  +  P3. 

In  the  two-repairmen  case,  the  repairmen  are  busy  when  the 
system  is  in  state  3: 
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For  P  ,   the  only  time  a  machine  must  wait  in  cold  standby 
is  when  the  system  is  in  state  0;  therefore,  P   =  p,. 

The  availibility  of  at  least  one  helicopter  in  the 
system  is  the  sum  over  all  states  n  of  the  percentage  of 
the  flyable  helicopters  in  each  state  multiplied  by  the 
proportion  of  time  the  system  is  in  that  state; 


Av  =  Pq(D  +  Px(2/3)  +  p2(l/3)  +  P3(0) 


The  estimate  of  the  availabiliy  is  calculated  by  substitut- 
ing p  ,  n  =  0,1,2,3,  into  the  above  expression. 
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III.   THE  SIMULATION 


A.   THE  COMPUTER 

The  microcomputer  used  for  this  simulation  was  the  IBM 
Personal  Computer.  Since  the  program  MACHREPR  uses  color 
graphics  extensively,  the  computer  requires  a  color  display 
monitor  and  either  IBM's  Color/Graphics  Adapter  or  a  suit- 
able graphics  card  from  another  manufacturer.  The  program 
will  not  run  without  a  graphics  adapter.  Other  minimum, 
hardware  requirements  for  this  simulation  are  128  kilobytes 
of  read-write  memory  (RAM)  and  at  least  one  54-inch,  floppy 
disk-drive . 

At  the  beginning  of  this  project,  the  author  chose  to 
use  IBM's  disk  operating  system  PC-DOS  with  the  intent  of 
enabling  the  simulation  to  run  on  other  microcomputers 
using  the  more  generic  version  of  PC-DOS:  MS-DOS  by 
Microsoft.  It  soon  became  apparent,  however,  that  MACHREPR 
would  need  to  directly  access  certain  parts  of  the 
computer's  memory  in  order  to  improve  the  program's 
performance.  Since  the  program  makes  direct  calls  to  the 
PC-DOS,  there  can  be  no  assurance  that  the  simulation  will 
run  on  other  "IBM-compatible"  machines. 


20 


B.  THE  PROGRAMMING  LANGUAGES 

MACHREPR  is  written  in  Advanced  BASIC  which  is  the  pro- 
gramming language  supplied  with  PC-DOS.  Advanced  BASIC  was 
selected  because  of  its  convenient  graphics  commands  and 
its  ready  availability.  Although  Advanced  BASIC  is  an 
interpreter,  program  speed  was  not  an  essential  criterion 
that  warranted  using  a  compiled  language.  Advanced  BASIC 
was  insufficient,  however,  in  two  aspects,  and  so  MACHREPR 
uses  two  subroutines  written  in  8088  assembly  language  to 
improve  program  performance.  These  two  subroutines,  named 
SCRNSHFT  and  RNGEM ,  will  be  described  later  in  this 
section . 

C.  PROGRAM  DESCRIPTION 

Appendix  B  contains  the  flowcharts  for  the  BASIC  pro- 
gram MACHREPR.  The  flowcharts  illustrate  the  simulation's 
general  construction  and  program  flow;  however,  they  are 
not  a  verbatim  duplication  of  the  program  listing.  The 
program  listing  for  MACHREPR  and  the  two  source-code  list- 
ings for  RNGEN  and  SCRNSHFT  are  in  Appendixes  C,  D,  and  E 
respectively . 

After  initializing  several  variables  and  arrays,  the 
program  presents  the  title  screen.  The  program  menu 
appears  next  with  the  following  options:  (1)  see  the  pro- 
gram instructions,  (2)  change  the  model's  parameters,  (3) 
use   the   default   parameters,   (4)   set   the   random-number 
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generator  seed,  or  (5)  end  the  program.  The  instructions 
that  are  available  on  the  screen  are  a  summarized  version 
of  those  found  in  Appendix  A,  the  user's  guide.  The  model 
parameters  that  can  be  changed  by  the  user  are  the  mean 
time  to  failure,  the  mean  time  for  repair,  and  the  number 
of  repairmen. 
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Figure  3:  An  Example  of  the  Graphics  Display. 

Figure  3  shows  an  example  of  the  graphics  display 
screen  as  it  appears  during  the  simulation.  The  upper  half 
of  the  screen  contains  the  system  state-versus-time  graph. 
The  system  state  values  are  on  the  ordinate  of  the  graph, 
and  fifty  divisions  representing  time  are  on  the  abscissa. 
The  three  graphic  figures  of  helicopters  in  the  lower  left 
corner  symbolize  the  status  of  the  system.   The  helicopter 
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seen  above  the  others  is  in  active  service;  the  two 
machines  with  broken  rotors  are  down.  In  this  example,  the 
helicopter  graphics  indicate  that  the  system  is  in  state  2, 
where  helicopter  B  is  flying,  helicopter  C  is  under  repair, 
and  helicopter  A  is  awaiting  repair.  The  graph  shows  that 
the  current  time  is  69  and  that  the  system  jumped  from' 
state  1  to  2  at  t  =  68. 

The  table  in  the  lower  right  corner  contains  the  esti- 
mated and  the  limiting  values  of  the  quantities  discussed 
in  Section  II.   The  row  labels  are: 

Ar   -  the   arrival   rate  of  machines   into  the   repair 
queue , 

Wq   -  the  average  wait  time  of  a  machine  in  the  repair 
queue , 

Lq   -  the   average   number   of   machines   in   the   repair 
queue , 

W   -  the  average  amount  of  time  a  machine  is  down, 

L   -  the   average  number  of  machines   in   the  repair 
system, 

Avl  -  the  proportion  of  time  that  at  least  one  machine 
is  available. 

Pd   -  the  probability  that  a  down  machine  must  wait  for 
repair , 

Ps   -  the  probability  that  an  up  machine  must  go  to 
cold  standby. 

MACHREPR  consists  of  four  major  parts:  (1)  main  pro- 
gram, (2)  clock  module,  (3)  failure  module,  and  (4)  repair 
module.  The  main  program  handles  the  housekeeping  details 
such   as   initializing   variables,   loading   the   machine 
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language  subroutines,  generating  the  graphics  display,  and 
performing  the  tasks  selected  from  the  program  menu.  Within 
the  main  program  are  routines  for  computing  the  limiting 
values  of  the  tabulated  quantities  and  for  printing  them  on 
the  display  screen.  The  main  program  also  calculates  the 
initial  failure  times  for  the  first  entry  into  the  clock 
module . 

The  clock  module  determines  the  next-event  time, 
whether  the  event  is  a  failure  or  the  completion  of  repair, 
and  which  helicopter  it  affects.  The  simulation  clock  is  a 
fixed-increment,  time-advance  type.  In  contrast  to  a  next- 
event,  time-advance  clock  which  skips  over  the  intervening 
time  between  events,  the  clock  module  in  MACHREPR  advances 
the  time  in  equal  increments.  With  this  type  of  clock,  the 
user  experiences  the  illusion  of  real-time  passage. 

During  the  execution  of  the  simulation,  the  clock 
module  monitors  the  keyboard  for  certain  commands  issued  by 
the  user.  Pressing  the  "P"  key  pauses  the  program  to  allow 
the  user  to  study  the  graphic  display.  Pressing  the  "C" 
key  causes  the  program  to  continue.  The  "S"  key  will  stop 
the  simulation  and  return  the  program  to  the  program  menu. 
Should  the  clock  reach  the  preset  limit  of  9950  time  units, 
the  program  will  automatically  stop  the  simulation  and 
return  to  the  program  menu. 

Every  time  the  clock  is  advanced  one  increment,  the 
clock  module   updates   the   system  state-versus-time  graph 


24 


according  to  the  current  state  and  time  values.  As  the 
time  advances,  the  plotted  lines  will  eventually  reach  the 
edge  of  the  graph.  When  this  condition  is  met,  the  clock 
module  shifts  the  entire  graph  left  one  time  unit  and  draws 
the  new  time  increment.  This  procedure  continues  for  all 
successive  time  increments.  When  the  simulation  clock  has 
reached  the  time  for  the  next  event,  the  program  jumps  to 
either  the  failure  module  or  the  repair  module  as 
appropriate . 

The  failure  module  changes  the  helicopter  display 
according  to  the  current  state  and  model  parameters .  The 
program  keeps  track  of  the  machines  individually,  so  that 
the  graphic  image  of  the  particular,  failing  helicopter 
transforms  from  a  green,  flying  figure  to  a  red,  grounded 
one.  If  a  repairman  is  free,  the  graphic  figure  of  a  man 
is  drawn  next  to  the  aircraft.  If  no  repairman  is  avail- 
able, then  the  human  figure  is  omitted,  signifying  that  the 
helicopter  is  placed  in  the  repair  queue.  If  a  repairman 
was  assigned,  the  failure  module  calls  the  random-number 
generator  and  computes  the  service  time.  The  sum  of  the 
service  time  and  the  current  time  is  the  time  that  the 
repair  will  be  completed  and  this  is  compared  to  the  next- 
event  times  of  the  other  machines  after  the  program  returns 
to  the  clock  module. 

In  a  similiar  manner,  the  repair  module  changes  the 
graphic  figure  of  the  helicopter  that  has  just  completed 
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repairs.  If  there  are  less  than  two  helicopters  flying 
when  the  repair  event  occurs,  then  the  repaired  machine 
transforms  from  red  colored  and  disabled  to  green  colored 
and  flying.  If  two  helicopters  are  already  flying,  then 
the  repaired  machine  goes  into  the  standby  queue.  The  now 
free  repairman  moves  to  the  helicopter  at  the  head  of  the 
queue,  or  if  the  queue  is  empty,  he  disappears  from  the 
screen . 

Before  the  program  returns  to  the  clock  module,  the 
repair  module  computes  the  interarrival  time  to  failure  of 
the  helicopter  just  repaired.  Failure  times  are  not  com- 
puted for  helicopters  placed  in  the  standby  queue.  If  the 
repairman  is  assigned  to  another  helicopter,  the  repair 
module  also  computes  the  service  time. 

After  returning  to  the  clock  module,  the  simulation 
updates  the  estimates  of  the  tabulated  quantities  and 
prints  them  on  the  display  screen.  Both  the  failure  and 
repair  modules  change  the  state-indicator  variable  to  its 
new  value  so  that  the  clock  module  will  draw  a  horizontal 
line  one  time  unit  in  length  at  the  proper  place  on  the 
graph . 

SCRNSHFT  is  the  subroutine  that  shifts  the  system 
state-versus-time  graph  one  time  unit  to  the  left  per  clock 
cycle.  The  subroutine  uses  certain  program  transfers 
called  interrupts  to  the  PC-DOS  Basic  Input/Output  System 
(BIOS)   which   contains   routines   that   control   the   video 
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display.  SCRNSHFT  is  not  appreciably  faster  than  a  BASIC 
routine  designed  to  perform  the  same  task;  however,  the 
subroutine  eliminates  an  annoying  flicker  in  the  display 
that  was  characteristic  of  the  BASIC  version. 

D.   PROGRAMMING  CONSIDERATIONS 

The  program  computes  the  interarrival  (failure)  and 
service  times,  assumed  to  be  continuous,  exponentially  dis- 
tributed, random  variables,  in  the  single-precision  format. 
The  current  time,  an  integer  value,  is  added  to  the  inter- 
arrival time  and  to  the  service  time  to  yield  the  next- 
failure  time  and  the  next-repair  time  respectively.  In  the 
unlikely  event  that  the  two  single-precision,  random 
variables  match,  the  simulation  will  proceed  as  if  the 
failure  occurred  earlier. 

Whenever  the  model  jumps  from  one  state  to  another,  the 
program  computes  a  next-failure  time  and  a  next-repair  time 
for  the  helicopter  just  affected  by  the  transition.  Thus 
MACHREPR  always  keeps  in  memory  six  next-event  times  for 
all  machines  regardless  of  their  individual  status  .  To 
determine  the  next  event,  the  clock  module  compares  the 
smallest  of  the  three  failure  times  with  the  smallest  of 
the  three  repair  times.  With  this  method,  there  must  be  a 
means  of  preventing  the  program  from  simulating  events  con- 
sidered impossible  by  the  assumptions  of  the  model. 
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For  example,  one  model  assumption  is  that  a  helicopter 
in  cold  standby  is  not  subject  to  failure.  It  follows  that 
the  same  machine  cannot  be  repaired  since  it  is  not  in  the 
repair  queue.  Similarly,  a  helicopter  that  is  in  repair 
cannot  fail  and  a  helicopter  that  is  flying  cannot  be 
repaired.  To  handle  these  impossible  events,  the  program 
will  set  the  failure  time  or  the  repair  time,  as  appropri- 
ate, to  equal  10    which  represents  infinity. 

Because  of  the  discrete  nature  of  the  graphic  display, 
all  time  values  must  be  converted  to  integers  prior  to 
plotting  the  state-versus-time  graph.  It  is  for  this  rea- 
son that  the  random  number  representing  the  next-event 
time,  be  it  a  failure  or  a  repair,  is  rounded  up  to  the 
next  higher  integer.  The  result  is  that  the  interarrival 
and  service  times  are  not  exponentially  distributed  but  are 
instead  geometrically  distributed. 

The  discreteness  of  the  time  variable  introduces  a  bias 
in  the  estimators  since  they  are  based  on  the  assumption  of 
continuous  distributions.  The  effects  of  the  bias,  how- 
ever, are  minimal.  The  values  of  the  estimates  are  still 
seen  to  converge  to  the  limiting  values  as  time  advances. 
Despite  the  generation  of  geometric  random  variables,  the 
simulation  still  provides  a  sufficient  visual  demonstration 
of  the  machine-repairman  model. 
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E.   RANDOM  NUMBER  GENERATION 

The  author  found  BASIC'S  random-number  generator,  a 
function  called  RND,  unsuitable  for  simulation.  Even  so,  a 
sample  of  10,000  uniform  variates  produced  by  RND  passed 
with  .95  confidence  three  standard,  nonparametr ic  tests  for 
uniformity  and  independence:  the  serial  test,  the  frequency 
test,  and  the  runs  up-and-down  test.  On  the  other  hand,  a 
two-dimensional  scatter  plot  of  5000  pairs  of  uniform  vari- 
ates, shown  in  Figure  4,  clearly  illustrate  the  lattice 
structure  of  the  RND  generator. 
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Figure  4:  Two-Dimensional  Plot  of  5000  Pairs  of  Uniform 
Variates  from  the  BASIC  RND  Function. 


MACHREPR  uses  a  random-number  generator  called  RNGEN . 
Written  in  8088  assembly  language,  RNGEN  is  modification  of 
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a  generator  programmed  by  Associate  Professor  B.O.  Shubert 
of  the  Naval  Postgraduate  School  for  use  on  the  IBM  PC. 
RNGEN  uses  the  algorithm  described  by  Fishman  [Ref.  4] 
which  is  based  on  the  expression 

X.^,  =  16807*X.  mod  (231  -  1) 
l+l  l 

where   X.   is   the   old   seed   and  X.  ,,   is   the   new  seed, 
l  l+l 

Dividing  the  new  seed  by  the  modulus  produces  the  uniform 
(0,1),  random  number. 

This  algorithm  is  the  one  used  by  LLRANDOM,  a  random- 
number  generator  developed  at  the  Naval  Postgraduate  School 
and  described  by  Lewis,  Goodman,  and  Miller  [Ref.  5]  as 
able  to  produce  good  quality,  uniform  variates . 

The  simulation  computes  exponential  random  numbers  by 
the  inverse  probability  integral  transformation  of  uniform 
variates  produced  by  RNGEN.  Let  U  be  a  uniformly  distri- 
buted, random  number  and  let  X  be  exponentially  distributed 
with  parameter  a.  An  exponential  random  number  is 
generated  by  setting  u  =  F(x)  and  solving  for  x, 

■at     ^         i     -ax 
u  =  F(xj  =1  -  e 

,         -ax 
1  -  u  =  e 

ln(l-u)  =  -ax 

x  =  -ln(l-u)/a,  or 

x  =  -ln(u)/a. 
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A.  INTRODUCTION 

MACHREPR  is  a  stochastic,  discrete-event  simulation  of 
a  machine-repairman  model,  which  is  a  special  case  of 
birth-death  processes.  The  model  consists  of  three  heli- 
copters, of  which  two  are  in  service  and  one  is  in  cold 
standby,  with  an  option  of  one  or  two  repairmen.  The  pro- 
gram output  is  a  graphics  display  composed  of  a  system 
state-versus-time  graph,  a  table  of  statistics,  and  animat- 
ed figures  that  illustrate  the  dynamics  of  the  process. 

This  manual  discusses  the  model  and  guides  the  user  in 
the  step-by-step  operation  of  the  program.  In  describing 
commands  and  keyboard  entries,  the  following  notation  is 
used.  Single  quotes  ('  ')  enclose  commands  and  phrases 
that  are  to  be  entered  verbatim.  The  act  of  pressing  a 
specific  key  is  indicated  by  brackets  (<  >)  encasing  the 
letter  or  letters  on  the  key.  For  example,  'format  b:  ' 
<ENTER>  means  type  the  command  exactly  as  it  appears 
between  the  quotes  and  follow  it  by  pressing  the  enter  key. 
This  user's  guide  shows  the  commands  in  the  lower  case; 
however,  the  IBM  PC  will  accept  upper  case  letters  as  well. 

B.  MODEL  DESCRIPTION 

MACHREPR  contains  two  versions  of  the  machine-repairman 
model:  the  one-repairman  case  and  the  two-repairman  case. 
In  both  cases,  the  nominal  mission  of  the  model  is  to  keep 
two   helicopters   flying   at   all   times.   If   the   third 
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helicopter  is  serviceable,  it  is  kept  in  cold  standby  as  an 
immediate  replacement  for  the  other  two.  Upon  failure  of  a 
helicopter,  a  repairman  is  assigned  if  he  is  available; 
otherwise  the  helicopter  joins  a  queue  to  await  repair. 

Disabled  (down)  helicopters  remain  in  the  repair  queue 
until  a  repairman  becomes  free.  The  repair  queue  is  a 
first  in,  first  out  type  (FIFO).  A  newly  repaired  helicop- 
ter will  go  immediately  into  flying  operation  if  less  than 
two  aircraft  are  in  service.  If  there  are  already  two 
helicopters  flying  when  a  repair  is  completed,  then  the 
third  goes  into  the  cold-standy  queue. 

There  are  several  simplifying  assumptions  for  these 
models.  All  helicopters  are  assumed  to  have  exponentially 
distributed  service  lives  with  the  same  failure  rate  a.  If 
T  is  defined  as  the  length  of  a  helicopter's  operating 
life,  then  the  mean  time  to  failure  (MTTF)  is  E(T)  =  1/or. 
The  simulation  uses  the  MTTF  as  the  input  parameter  that 
establishes  the  particular  exponential  distribution  for  the 
operating  life.  Helicopters  will  fly  continuously  until 
failure.  A  machine  that  is  in  cold  standby  is  assumed  not 
to  be  subject  to  failure. 

Another  assumption  is  that  the  service  time  for  the 
repairmen  has  an  exponential  distribution  with  parameter  y. 
If  S  is  the  time  required  to  complete  a  repair  job,  then 
the  mean  time  for  repair  is  E(S)  =  1/y.  Like  the  MTTF,  the 
mean  time  for  repair  (MTFR)  is  an  input  variable  for  the 
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simulation.   The  repairmen  are  assumed  to  be  equally  compe- 
tent and  thus  have  the  same  MTFR. 

In  both  cases  of  the  machine-repairman  model,  the  state 
of  the  system  is  defined  as  the  number  of  helicopters  down 
at  time  t.   The  state  variable  is  N(t)  =  n,  n  =  0,1,2,3. 

C.   HARDWARE  AND  SOFTWARE 

1 .  Hardware 

The  simulation  is  programmed  for  the  IBM  Personal 
Computer  (IBM  PC)  with  a  color  monitor.  Since  the  program 
uses  color  graphics,  the  microcomputer  must  be  equipped 
with  IBM's  Color/Graphics  Adapter  or  a  suitable  graphics 
card  from  another  manufacturer.  MACHREPR  will  not  run 
without  the  graphics  adapter.  The  program  uses  two 
assembly-language  subroutines  that  are  loaded  in  memory 
outside  of  BASIC'S  64  kilobyte  (K)  workspace;  therefore, 
there  must  be  a  minimum  of  128  K  of  memory  installed  in  the 
microcomputer.  The  computer  must  have  at  least  one 
54-inch,  floppy  disk-drive. 

2 .  Software 

MACHREPR  and  the  two  subroutines  are  supplied  on  a 
distribution  diskette.  The  user  must  provide  the  disk 
operating  system  PC-DOS  2.0  (or  2.10)  and  Advanced  BASIC 
(BASICA)  in  order  to  run  the  simulation.  MACHREPR  uses 
BASIC'S  POKE  command  and  certain  program  transfers,  called 
interrupts,  to  the  Basic  Input/Output  System  (BIOS)  in  the 
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IBM  ROM  (read-only  memory).   Because  of  these  interrupts, 
there  can  be  no  assurance  that  the  simulation  will  operate 
properly  on  other  "IBM  compatable"  machines. 
3 .   Program  Files 

The   distribution   diskette   contains   the  following 
program  files: 

*  MACHREPR.BAS  -  the  main  program. 

*  RNGEN.SRT  -  a  random-number  generator  written  in  8088 
assembly-language  used  as  a  subroutine  to  the  main 
program. 

*  SCRNSHFT.SRT  -  another  assembly-language  subroutine 
used  to  perform  a  specific  task  on  the  graphics 
screen . 

*  RNGEN.LST  -  a  source-code  listing. 

*  SCRNSHFT.LST  -  a  source-code  listing. 

*  MACHREPR.BAT  -  a  batch  file  that  loads  and  starts  the 
simulation . 

The  simulation  requires  the  first  three  files  listed  above 

in  order  to  run.  The  two  subroutines  are  binary-image  files 

and  can  not  be  executed  from  PC-DOS.   The  two  source-code 

listings   are   provided   as   a   convenience   to   the   user. 

MACHREPR.BAT  is  a  simple  example  of  how  a  batch  file  can  be 

used  to  easily  start  the  simulation. 

D.   GETTING  STARTED 

1 .   Making  a  Backup  Copy 

The  first  step  is  to  make  a  copy  of  the  distribu- 
tion  diskette.   After   the   copy   is   made,   the   original 
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diskette  should  be  stored  in  a  safe  place  and  used  only  to 
make  a  replacement  for  an  unserviceable  application 
diskette.  Although  there  are  several  ways  that  the  user  may 
choose  to  create  an  application  diskette,  this  manual  will 
suggest  just  two: 

*  A  self -starting  (bootable)  diskette  containing  the 
operating  system  files,  the  BASICA  command  file,  and 
the  simulation's  program  files. 

*  A  simple  backup  copy  which  contains  only  the  program 
files  and  which  must  be  used  in  conjunction  with  a 
system  diskette. 

Owners  of  microcomputers  with  a  single  disk-drive  may  find 

the  first  option  more  convenient. 

To  make  a  bootable  diskette,  follow  the  instruc- 
tions in  the  PC-DOS  manual  and  format  a  blank  diskette 
using  the  FORMAT  command  with  the  /S  parameter.  The  /S 
parameter  causes  the  computer  to  transfer  the  operating 
system  files  from  the  DOS  diskette  to  the  newly  formatted 
diskette.  The  Advanced  BASIC  file,  called  BASICA.COM,  is 
also  found  on  the  DOS  diskette  and  can  be  copied  onto  the 
application  diskette  with  the  COPY  command.  Finally,  use 
the  COPY  command  again  to  transfer  the  simulation's  six 
program  files  from  the  distribution  diskette  onto  the 
bootable  diskette. 

If  the  second  option  is  desired,  then  format  a 
blank  diskette  without  the  /S  parameter.  This  procedure 
does  not  transfer  the  operating  system  files  to  the  appli- 
cation  diskette.   The   next   step   is  to  use  the  DISKCOPY 
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command   to   make   an   identical   copy   of   the   distribution 
diskette . 

2  .   Starting  the  Simulation 

If  the  first  option  was  used,  insert  the  bootable 
diskette  into  drive  A  and  turn  on  the  microcomputer.  After 
a  few  seconds,  the  DOS  prompt  A>  will  appear  on  the  screen. 
It  is  here  that  the  batch  file  MACHREPR.BAT  can  be  used. 
Simply  type  'machrepr'  <ENTER>  to  start  the  simulation. 

If  using  the  non-bootable  diskette,  then  turn  on 
the  microcomputer  with  the  DOS  diskette  in  drive  A.  Next 
load  Advanced  BASIC  into  memory  by  typing  'basica'  <ENTER>. 
After  the  computer  is  in  BASIC'S  command  mode,  remove  the 
DOS  diskette  from  drive  A  and  insert  the  diskette  contain- 
ing the  simulation.  Type  'load  "machrepr . bas  ",  r  '  <ENTER> 
to  bring  the  program  into  memory  and  to  start  the 
simulation  . 

If  desired,  the  simulation  can  be  started  from  a 
second  disk-drive.  Begin  as  before  by  turning  on  the 
microcomputer  with  the  DOS  diskette  in  drive  A  and  loading 
BASICA.  With  the  application  diskette  inserted  into  drive 
B,  type  'load  "b rmachrepr . bas  "  , r '  <ENTER>.  During  execu- 
tion, MACHREPR  reads  the  program  diskette  and  expects  to 
find  it  in  the  default  disk-drive.  When  the  program  fails 
to  find  the  required  files  on  the  DOS  diskette  in  drive  A, 
it  asks  the  operator  for  the  correct  drive  label.  Simply 
respond  to  the  program's  question  by  pressing  <b> . 
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Figure  5:  The  Title  Screen. 

3.   Title  Screen 

The  display  will  remain  dark  for  a  few  seconds 
while  the  program  initializes  variables,  arrays,  and  func- 
tion definitions.  When  initialization  is  done,  MACHREPR 
will  present  the  title  screen  as  shown  in  Figure  5. 

E.   THE  PROGRAM  MENU 

The  program  menu,  depicted  in  Figure  6,  appears  next 
and  provides  the  means  of  controlling  the  simulation.  With 
the  options  available  on  the  menu,  the  user  may  choose  to 
read  the  program  instructions,  change  the  model  parameters 
or  use  the  default  model  parameters.  The  program  instruc- 
tions that  are  available  on  the  screen  are  a  summarized 
version   of   those   found   in   this   manual.   The   model 
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parameters  that  can  be  altered  by  the  user  are  the 
helicopter's  mean  time  to  failure,  the  repairmen's  mean 
time  for  repair,  and  the  number  of  repairmen  to  employ. 
Additionally,  the  initial  seed  for  the  random-number 
generator  can  be  set. 
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Figure  6:  The  Program  Menu. 

The  simulation  starts  and  ends  with  the  program  menu. 
Whenever  the  user  stops  the  simulation  or  whenever  the  sim- 
ulation clock  reaches  the  preset  limit  of  9950  time  units, 
the  program  returns  to  the  menu.  The  simulation  then  can 
be  restarted  with  perhaps  a  difference  set  of  model  parame- 
ters, or  the  program  can  be  ended  altogether. 
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1 .  Setting  the  Random-Number  Generator  Seed 

To  set  the  seed,  press  the  <s>  key.  The  program 
will  then  respond  by  asking  for  the  value  of  the  seed.  The 
seed  is  entered  by  typing  any  integer  in  the  range  of  1  to 
2,147,483,646  (=  231  -  2)  and  pressing  <ENTER>.  The  pro- 
gram will  accept  a  non-integer  value;  however,  during 
processing  the  number  will  be  truncated  to  an  integer,  with 
the  result  possibly  being  different  than  the  desired  value. 
After  the  seed  has  been  entered,  the  program  returns  to  the 
program  menu. 

It  is  not  necessary  to  set  the  random-number  gener- 
ator seed  to  run  the  simulation.  There  is  a  default  seed 
embedded  in  the  generator's  program  code,  and  it  is  updated 
automatically  with  each  call  for  a  random  number.  If  the 
simulation  is  stopped  and  then  restarted  without  setting 
the  seed,  the  generator  will  use  as  the  initial  seed  the 
value  remaining  from  the  last  call  in  the  last  simulation 
run . 

2.  Changing  the  Model  Parameters 

When  <c>  is  pressed,  MACHREPR  will  successively  ask 
for  the  desired  values  of  the  mean  time  to  failure,  mean 
time  for  repair,  and  the  number  of  repairmen.  The  quanti- 
ties are  entered  by  typing  a  positive  integer  and  pressing 
<ENTER>  in  response  to  each  question.  For  best  results, 
the  values  for  the  mean  time  to  failure  and  the  mean  time 
for  repair  should  be  at  least  ten.   The  program  will  accept 
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and  function  properly  with  non-integer  values  of  MTTF  and 
MTFR ;  however,  these  values  will  be  printed  on  the  graphics 
display  as  integers.  Negative  numbers  will  not  be  accepted 
for  the  MTTF  and  the  MTFR.  The  digits  1  and  2  are  the  only 
valid  responses  for  the  number  of  repairmen.  The  program 
will  immediately  start  the  simulation  after  the  number  of 
repairmen  has  been  entered. 

3  .   The  Default  Parameters 

If  this  option  is  selected,  then  the  simulation 
starts  with  the  model  parameters  set  at  default  values. 
These  values  are: 

MTTF  =  15  time  units, 

MTFR  =  10  time  units, 

Number  of  repairmen  =  1. 
The  default  parameters  are  selected  by  pressing  <d> . 

4.   Ending  the  Program 

To  end  the  program,  press  <e>.  MACHREPR  will  clear 
the  screen,  exit  the  graphics  mode,  and  leave  the  computer 
in  the  BASIC  command  mode.  If  the  user  wants  to  exit  BASIC 
and  return  to  the  DOS,  type  'system'  <ENTER> . 

F.   THE  GRAPHICS  DISPLAY 

MACHREPR  is  designed  to  give  the  observer  the  illusion 
of  real-time  passage  as  the  probabilistic  events  occur  in 
the  machine-repairman  system.  The  graphics  display  indi- 
cates the  model's  activities.   As  shown  in  Figure  7,  the 
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graphics  display  is  divided  into  three  sections:  the  state- 
versus-time  graph,  the  helicopter  display,  and  the  table  of 
quantities  used  with  queueing  models. 
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AR          0.132          0.092 
Wq          5.444         6.691 
Lq          0.721          0.617 
W          13.444       14.691 
L             1.485          1.355 
Avl        0.505          0.548 
Pd          0.529          0.458. 
Ps          0.235          0.262 

Figure  7:  An  Example  of  the  Graphics  Display. 

1.   The  State-Verses-Time  Graph 

The  upper  half  of  the  display  contains  the  state- 
versus-time  graph.  The  values  of  the  state  variable  N(t), 
the  number  of  machines  down  at  time  t,  are  on  the  ordinate. 
There  are  fifty  divisions  on  the  abscissa  which  represent 
units  of  time.  The  current  time  of  the  system  is  printed 
at  the  right-hand  end  of  the  abscissa. 

With  each  unit  advancement  of  the  simulation  clock, 
the  program  plots  a  horizontal  line  one  time  unit  in  length 
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at  the  height  corresponding  to  the  current  state  of  the 
system.  The  accumulation  of  these  individual  plots  forms  a 
step  graph  as  shown  by  the  example  in  Figure  7.  When  the 
graph  reaches  the  end  of  the  display,  the  program  makes 
room  for  the  next  plot  by  calling  SCRNSHFT  to  shift  the 
entire  graph  one  unit  to  the  left.  In  this  way,  the 
model's  activities  for  the  past  fifty  time  units  can  be 
seen  in  a  single  glance. 

2.   The  Helicopter  Display 

The  three  helicopter  images  in  the  lower  left 
corner  of  the  graphic  display  indicate  the  current  state  of 
the  system.  The  display  in  Figure  8  shows  how  the  helicop- 
ters are  arranged  in  the  initial  state  0.  All  three 
machines  are  colored  green.  Helicopters  B  and  C,  seen  in 
the  upper  half  of  the  display,  are  flying.  Helicopter  A  on 
the  "ground"  is  in  cold  standby. 

Since  all  helicopters  are  initially  either  flying 
or  in  cold  standby,  the  next  event  will  be  a  failure  in 
either  machine  B  or  C.  A  machine  that  fails  is  grounded 
and  transformed  into  a  red  colored  helicopter  with  a  broken 
rotor.  Because  there  is  only  one  down  helicopter,  the  sys- 
tem is  in  state  1  and  so  there  is  at  least  one  idle 
repairman.  The  program  draws  a  graphic  figure  of  a  man 
beside  the  down  aircraft  to  represent  a  machine  that  is  in 
repair.  For  example,  Figure  9  depicts  state  1  with 
helicopter  B  in  repair. 
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Figure  8:  Initial  Helicopter  Display  in  State  0 
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Figure  9:  Helicopter  Display  in  State  1. 

As  the  simulation  proceeds,  random  numbers  are  gen- 
erated and  transformed  into  interarrival   (failure)   times 


47 


and  service  (repair)  times.  The  program  changes  the 
appropriate  graphic  figures  when  the  simulation  clock 
reaches  the  time  for  the  next  event. 

Figure  10  shows  an  example  of  a  helicopter  display 
that  indicates  state  2.  In  this  case,  C  is  in  repair. 
Helicopter  A  is  also  disabled,  and  since  it  does  not  have  a 
repairman,  A  is  in  the  repair  queue. 
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Figure  10:  Helicopter  Display  in  State  2. 

In  Figure  11,  all  three  machines  are  down,  thus  the 
system  is  in  state  3.  There  is  a  repair  queue  of  length 
two  in  this  example. 

It  should  be  noted  that  Figures  10  and  11  are 
illustrations  from  the  one-repairman  model  since  only  one 
repairman  appears  in  either  display.   If  the  simulation  is 
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running  the  two-repairman  model,  then  the  program  will  draw 
the  second  repairman  beside  the  helicopter  that  fails  after 
the  one  already  in  repair. 
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Figure  11:  Helicopter  Display  in  State  3. 

3 .   The  Quantities  of  Interest 

In  the  study  of  queueing  systems,  there  are  several 
quantities  of  interest  that  are  used  to  evaluate  the  sys- 
tem's performance.  In  the  lower  right  corner  of  the 
graphics  display  (see  Figure  7)  there  is  a  table  of  the 
more  commonly  used  quantities. 

The  row  labels  have  the  following  definitions: 

*  AR  -  The  average  arrival  rate  of  machines  into  the 

repair  system. 

*  Wq  -  The  average  waiting  time  of  a  machine  in  the 

repair  queue. 
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*  Lq  -  The   average   number  of  machines   in  the  repair 

queue . 

*  W   -  The  average  amount  of  time  that  a  machine  is  in 

the  repair  system.   This  quantity  is  the  sum  of 
Wq  and  the  expected  service  time. 

*  L   -  The   average   number  of  machines   in  the  repair 

system.   This   quantity   includes  those  machines 
in   the   repair   queue   and   those   machines   in 
■repair . 

*  Avl-  The    machine    availability   defined    as    the 

proportion   of   the   time   that   at   least   one 
helicopter  is  flyable. 

*  Pd  -  The  probability  that  a  down  machine  must  wait 

for  repair. 

*  Ps  -  The  probability  that  an  up  machine  must  wait  in 

cold  standby. 

The  numbers  that  appear  in  the  right-hand  column 
under  the  heading  "Limit"  are  the  limiting  values  of  the 
quantities  of  interest.  In  other  words,  the  quantities 
converge  to  the  limiting  values  as  time  goes  to  infinity. 
The  limiting  values  are  calculated  with  the  balance  equa- 
tions that  correspond  with  the  particular  machine-repairman 
model . 

The  values  in  the  left-hand  column  are  estimates  of 
the  limiting  values.  These  quantities  are  computed  with 
data  collected  from  the  simulation.  Upon  the  occurance  of 
an  event,  the  program  recomputes  the  estimates  after  taking 
into  account  the  newly  collected  data  from  the  last  state. 

Because  the  plots  on  the  state-versus-time  graph 
must  be  given  in  terms  of  integer  coordinates,  all  time 
variables  must  be  rounded  up  to  the  next  higher  integer.  As 
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a  result,  the  estimators  are  biased  since  they  are  based  on 
the  assumption  of  continuous  time  when  in  fact  they  are 
calculated  with  data  based  on  discrete  increments  of  time. 
The  effects  of  the  bias  are  minimal;  the  simulation  still 
provides  a  sufficient  demonstration  of  the  convergence  of 
the  estimates  to  the  limiting  values. 

G.   KEYBOARD  COMMANDS 

There  are  three  keyboard  commands  available  to  the  user 
as  the  simulation  is  in  progress.  When  <p>  (for  pause)  is 
pressed,  the  simulation  freezes  to  give  the  user  the  oppor- 
tunity to  study  the  graphic  display.  To  continue  the  simu- 
lation, press  <c> .  The  third  command  is  <s>  for  stop  the 
simulation.  The  <s>  key  will  send  the  program  back  to  the 
program  menu. 

H.   PROGRAM  LIMITATIONS 

MACHREPR  always  starts  a  simulation  under  the  same  con- 
ditions. The  initial  state  is  N(0)  =  0  with  helicopters  B 
and  C  in  flight  and  helicopter  A  in  cold  standby.  There 
are  no  provisions  in  the  program  for  starting  the 
simulation  in  another  state.  Restarting  the  simulation 
from  the  program  menu  causes  the  current  time  and  all  data 
collection  variables  to  be  set  to  zero. 

As  noted  before,  all  event  times  are  rounded  up  to  the 
next   higher   integer   for   plotting   the   state-versus-time 
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graph.  With  exponentially  distributed  interarrival  and 
service  times,  there  is  a  non-zero  probability  that  the 
sojourn  in  any  particular  state  will  be  less  than  one  time 
unit.  The  program  will  still  alter  the  helicopter  display 
properly  and  compute  the  estimates  correctly;  however,  the 
state-versus-time  graph  may  display  a  jump  of  two  states 
because  of  the  program's  inability  to  plot  a  line  of  less 
than  one  time  unit  in  length.  The  number  of  times  this 
discrepancy  occurs  can  be  minimized  if  the  MTTF  and  the 
MTFR  are  set  at  values  greater  than  ten. 
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APPENDIX  B 
MACHREPR.BAS  FLOWCHART 
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QBSUI   1021Q 


Ftnd  *htcn  helo 

•111   fall   next: 

FARINGHELOl 


GQSUfl  UJ22H 


Ftnd  *ntcn  nelo 

f*  reoatred: 

REPAIREDHELOS 


Clock   Module 


Set  nelo  A  fall 
time  to  Infinity 


Commute  B's 
fafiwe  tine 


Co«out«  C*» 
fat lur«  ttne 


S«t  all  repair 

tines  to  Infinity 


STOSYHELOi  -  'A* 


/ WHILE 
-*<f  NXTEVNT<995 


STATE-0  ? 


STATE-3  ? 


Set  NXTFAIL 
Infinity 


NXTREPR  - 
MIN(REPRTIM£(.)) 


GOSUB  1Q??Q 


Ftnd  wntcn  nelo 

1«  reDafred: 

REPAIREOHELOt 


NXTEVNT  »  NXTREPR 


NXTFAIL    - 
KIN(FAILTIME( 


)) 


Set  NXTREPR 
Infinity 


QPSUB  1Q21Q 


Find  ■ntcn  nelo 

»111  fafi  next: 

FAILINGHELOj 


NXTEVNT    -   NXTFAIL 
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WHILE 
T1WE<NXTEVNT 


uooate  TlM€ 


CM.L  SCRNSHFT 


Shfft  U«« 

gT-«DM  to  left 

one  unit 


Check  keyOoard 
(Xrffer 


Pauee  routine 


Add  sojourn   ttn«    fn  stats   to 
cvawui a t f v»   U»e    In  state 


LASTEVNT-NxTEVNT 


Gosua  ™aa 


REPAIR 
Module 


CP$UB  2QQQ 


FAILURE 
Module 


uooate 

statistics 

and  D«"tnt   to 

screen 


Cl 
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Clock  Module  Suorouttnes 


rSTwrr 
SUroutlr 


fAIUHOClOt  • 


F4IUNQC10*  •   I 


FAIUMOCLOS  - 


EePAnraHoot 


REP*nSOH€LOt 


/"     5T/U»r      ^\ 
/     SUirouttn*     A 

" 

I    "™    ) 

/      pqixt          / 

' 

,  ] 

! 

/  Qnm  rut»-      / 
/     «rw  tt—     / 

/          flrW         / 

! 

Star*  n*io  A   In 
•my  OfiMCLO 

' 

J 

' 

' 

fOrtm  StAlC  -  •/ 
/      rwlfcoettr*   / 

Star*  antlr* 

wwn   In  «rr»y 

MAfMS09( 

/                             / 

..       } 

f 

/^RET 

urn"\ 
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FaMure  Module 


At   ■    3 

X   -    107 


[Subrout fne 
2000       J 


Iocre««nt 

SUHTAIL 


Incrnn»nt   STATE 


F2 


KS 


At   -    ? 


X    -    54 


cosua  1210Q 


Draw  4  «tor» 

r*d  r>«lo 

array* 


FIRSTrAILUPt-fALSE 


F3 


At   •    1 


X    -     1 


F2 


^/ 
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CALL    RNCEMU) 


REPRTlMe(AS)-TIMe 
••fNEXPONTCU.HU) 


rAILTIHe(AX)    -    JE31 


AX 


107 


cosua   12228 


Ground  1 lying 

it«lo    A    astlfjn 

repatfan 


CALL   RNC£N(U) 


FAlLTIME(At)-nME 
*FNEXPONT(U.ALP«A) 


Failure  Module  -  Part  A 


F3 

V 


At    -    1 


STDBYHELOt-null 


CJSUa    1??30 


Mov«  standby 

helo    to   1 lyfng 

statu* 


RETURN 
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casua   13?1Q 


Q-oona   flying 
h«  1 o :   no 
repairman 


fAIU71HC(AS)    -    1E31 


REPRTlME(At)    -    1E31 


WAITREPRt   -    *C 


RETURN 


Failure  Module  -  Part   B 


WAITREPBS    -    'A" 


IILINGHCLOS^X *JUAITREPRS   -    'B* 


man  1222a 


Grovnd  n«lo   A 
s««fgn  riMfrain 


CALL   RfKZMU) 


REPRTIHE(AI)-T1« 
*FHexPOM(U.HU) 


FAJLTIK(At)    -    1E31 


WAITREPRJ   -   null 


RETURN 
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F2 
c 


Failure  Module  -  Part  C 


UEPRTIMEtAI)    -    1E31 


■A}LTIMC(AI)    -    1E31 


MUMRMAN 
2? 


WAITR£PR$    -    *C 


UAITRCPflS   -    'A' 


WAI7REPR$   -    'B* 


CQSUB   122UL 


GrotryJ   flyfng 
n»lo:   no 
rep<  Iritan 


RETURN 
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Fat  lure  Module  Subroutines 


STMT 

SuOrowttn* 

12100 


STAflT 
SUrouttf 
13129 


CAU   flNCTN(U) 


A|    •    T1W   ♦ 


x  -  Mi  il  •  it 


RTPOTIKC})   -  At 


r»iux«(j)-mi 


I  ■   lift   XI   ■    131 


OfWTIK(J)  -  Al 


FAUTrxiu-icJi 


Wd.115).   GSMC1.0 


Otm  r«d  n«io 
•nd  «tor«  In 
•rrty  fltOHCVO 


Ortm  rtottraan 
•nd  rtort  «»th 

h»'o  to  »rr»y 


RETURN 


pum.ios).g>w€io| 


puT(x,i4>),acPK\o 


0«'«y  routfrwj 

coum   Irvm  1   to 

3M 


PUT(».lA3)tBTDMCL0 


PUT(>t143)tarPHO,M 


RETURN 
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Repafr  Module 


[Subroutf ne 
3000       J 


AI  •   3 

X   -    107 


At   -    2 


X    -    54 


-    2 


Rb 


At   -    1 


t    -,    1 


Rc 
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STD«Y>CLOf  -   •*• 


Repair  Module  -  Part  A 


ON  AS 

GOTO 


noiYHcio*  -  •!• 


<> 


FAllTIIC(AI)   -   1E31 


WPflTDC(Al)   •    1CJ1 


com  nan. 


Owngt  rod  holo 
to  cold  stanaoy 


STATE  -  STATE  -   1 


$TO»T>CLOf  •   •C» 


RETURN 
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Repafr  Module  -   Part   B 


CALL  RNCEMU) 


M]UJ«(AI)   -  1M 
'  *«>POMT(U. ALPHA) 


,.t 


FP«TIJ<{AI)    .    mi 


arc*  '*"■ 


ffvang*  htlo  ta 
groan  and  fly 


MUMWAN  •   1 


state  •  STATE  •  1 


RETURhT\ 


J 


AS   -    )S    X 


»I   •    it    »   •    ) 


<y 


Add  waft  tla* 

in  quoua  to 

total   of  aait 

tfaw*  In  quaua. 


GQSUI   13038 


CTianga  rtd  halo 
to  ran  hwlo 


CALL   RNCCM(U) 


RCPRTI)C(Af)   •   TIME 
♦   PKXPONT(Utm) 


WAITHtPQl   •  •• 


STATf    .    STATE    -    1 


RETURN 
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Repafr  Module  -  Part  C 


CALL   RWGEN(U) 


FAIT.TIX{it)    -    TIME 
*    F*€XPONT(U, ALPHAS 


*EP«TD*(AJ)    -    1E31 


sraia   i3o?o 


Cr>«ng»   to   X 
hmto  tnd  f I 


»y 


AiS  -  l:  x  -  l 


AiX   -   2:    X   -    5* 


A1J  -   3:    X   -   107 


A<H  Mft   tt«» 
?n  Qucut   to 
total   weft  tta* 
tn 


co^ia    13010 


cnang#  r*d  n*]o 
to  r*c  n«lo 


CALL   RWGEX(U) 


R£P€PfiTI>«(AlX)   - 

TIHE    ;    FS£XP<XT(U.HU) 


i[ 

VATTRfPB*    -   wm 

" 

STATF    -   STATE    -    1 

VAITREPtfS  -  'B' 


STATE    -    STATE   -    1 


RETURN 
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Repair  Module  Subroutines 


STMT 

Sobroutln* 

13010 


PUTQt.m)  BC0»€103 


I  wjT(»,io)t  cawcio 


RETURN 


START 

Subroutine 

1J0I0 


OoUy  reutfrwt 
count  froa  1  to 


200 


1  puTQt.io).  acoMCLO 


PUT(»,141),agTJKIM 


rSTMT 
Sufirouttn* 


13030 


PuT(X,lQ),gXTXlCT  ] 


PUTOt.lO),   CgjgCLg 


D«l«y  rowttntt 

count  froa  1  to 

MO 


PUT(X<1«3),    CBNCIO 


1  POT(X.IOS).   CWCL0~1 


RETURN 


RETURN 
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APPENDIX    C 
MACHREPR.BAS    PROGRAM    LISTING 


0  -MACHREPR.BAS,  Ver  3.1,  R.  E.  Nelsen,  30  Aug  1984 

1  CLEAR  :  KEY  OFF:   CLS 

2  GQSUB  30000:  'Check  for  BASICA  and  color/graphics  adaptor 

3  OEFINT  A-Z 

4  DIM  P'(3)1SUHP!(3),6RNHEL0(26B),REDHEL0(268),REDHEL02(268> 

5  DIM  TIMECHRT(1609),HAINSCRN<8002) 

6  DIM  FAILTIME'(3),REPRT1HE,(3),STARTQUE(4) 

7  RN6EN=0:TRUE=-1:FALSE=0:U'=0:SCRNSHIFT=3B4 

8  FMTl$='i#8.#it":FHT2$="tttfiDRVE$="A:' 

9  FIRSTSCRN=TRUE:WAtTREPR*="":FAILINGHELQS="":STDBYHELO$=" 

10  REPAIREDHEL0$=":S6$  =  STRING* (40,223) 

50  DEF  FNEXP0NT! (Al ! fA2! )  =  -  L0G(A1')/A2! 

51  DEF  FNROUNDUP(A')  =  INTIA!)  +  1 

52  DEF  FNHIN!(A!,B!,C')=(-(A!<BI  AND  A!<C!)«A!)+HB!<A!  AND  B!<C!)«B!>+(-(C!<A!  AND  C!<B!)*C!) 

53  DEF  FNHEL0P0SIT(A$)=-(  INSTRCA8C",A$)=l)  +  <-<  INSTRI"ABC",A$)=2)«54)+l-(  INSTR("ABC",A$)=3)M07) 
55  ' 

60  GOSUB  20000:  'Title  screen 

80  ON  ERROR  GOTO  30020:  'Load  randoi  ninber  generator. 

82  DEF  SEG  =  iHlAOO:  BLOAD  DRVE*+"RNGEN.SRT",0 

84  ON  ERROR  GOTO  0 

85  ' 

90  BLOAD  DRVE$+'SCRNSHFT.SRT",3B4:  'Load  screen  shift  subroutine 

91  ' 

100  '=====MAIN  PROGRAM  MENU===============— ================ 

120  SCREEN  0,1:C0L0R  14,3,0:MIDTH  40:CLS 

121  LOCATE  7,14:PRINT  "PROGRAM  HENU":PRINT  SB* 

122  LOCATE  9,1:PRINT  "(Instructions." 

123  LOCATE  11,1: PRINT  "(Dhange  aodel  parameters." 

124  PRINT:PRINT  "<D>efault  lodel  paraieters." 

125  PRINT:PRINT  "<S>et  the  randot  nueber  generator  seed." 

126  PRINT:PRINT  "<E>nd  the  program":  PRINT  SB$ 

130  LOCATE  21,2:  PRINT  'Enter  your  selection..."; 

131  ' 

140  GOSUB  10000:AI=  INSTR("ICDSEicdse",Af ) 

141  IF  AI=0  THEN  120  ELSE  ON  Al  6QT0  20100,200,210,150,900,20100,200,210,150,900 

142  ' 

149  '=====RAND0M  NUMBER  6ENERAT0R  SEED  R0UTINE============ 

150  COLOR  14,0,0:CLS 

151  LOCATE  1,3:  PRINT  "*  SET  RANDOM  NUMBER  6ENERAT0R  SEED  «":  PRINT  SGJ 

152  LOCATE  4,2:PRINT  "Permissible  seed  values  are  integers" 

153  PRINT  "  in  the  range:  1  to  2147483646." 

154  LOCATE  7,2:INPUT  'Enter  the  seed  value. ..";AI 

155  IF  Aid  OR  AD2147483646I  THEN  LOCATE  7,2:  PRINT  STRING* (39, 32> :  6OT0  154 
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159  DEF  SEG  =  &H1A00 

160  Al!=  INT  (Af/16777210I>:A2!=  INT ( CAI-A1 ' *16777210i) /i5536 ! ) 
142  'Poke  the  seed's  upper  2  bytes  into  RNGEN's  seed  storage. 

163  POKE  &H164,A1!:  POKE  M163,A2' 

164  A*=AI-Al!fl6777210i-A2!t65536!:Al'=  INT(At/256):A2!=  AI-AIH256 

165  'Poke  the  seed's  lower  2  bytes  intD  RNGEN's  seed  storage. 

166  POKE  iH162,Al!:  POKE  iH161,A2! 

167  GOTO  120:  'Return  to  ienu 

168  ' 

199  -=====PflRAHETERS  R0UTINE=============================== 

200  COLOR  15,5,0:CLS:HSG$=">>Enter  a  positive  integer  only.<<" 

202  LOCATE  1,6:PRINT  "*  CHANGE  MODEL  PARAMETERS  »":PRINT  SG* 

203  LOCATE  4tl: INPUT  'Enter  Mean  Tiie  To  Failure. ..",A!:IF  A'<=0  THEN  LOCATE  23,3:PRINT  MSG$;:LOCATE 
4,1:PRINT  STR1NG$(40,32);:G0T0  203  ELSE  ALPHA!=1/A' 

204  LOCATE  6,1: INPUT  "Enter  Mean  Tue  For  Repair...', A!: IF  A!<=0  THEN  LOCATE  23,3:PRINT  MS6$;:L0CATE 
6,1:PRINT  STR1N6*(40,32);:G0T0  204  ELSE  HU'=l/A! 

205  LOCATE  23,1:PRINT  STRING* <40 ,32) ;:HSG$=" 

206  LOCATE  8,1: INPUT  "Enter  nuaber  of  repainen  (1  or  2)...  ",NUMRMAN:IF  (NUMRHANOl)  AND  (NUMRHAN02 
)  THEN  LOCATE  B,1:PRINT  STRING*  (40 ,32) ;: GOTO  206  ELSE  250 

210  ALPHA!=l/t5:MU'=l/10:NUHRMAN=l 

250  NXTFAIL!=0:NXTREPR!=0:FIRSTFAILURE=TRUE:NXTEVNT=0:TIME=0:SUHFAIL=0:LASTEVNT=0:STATE=0:XX=31:SUHH 

Q=0 

252  FOR  1=0  TO  3:SUMP! ( I >  =0: STARTQUE { I ) =0: NEXT 

260  ' 

299  '=====THEORETICAL  STATISTICS  R0UTINE================== 

300  P!(1)=2*ALPHA!/MU! 

320  IF  NUMRMAN=1  THEN  P! (2)=P! (1)»P! (l):P! (3)=P! (2)*ALPHA!/MU!  ELSE  P! (2)=P! (1)*ALPHA!/HU!:P! (3)=P! ( 
2)*ALPHA!/(MU'+MU!) 

324  P!(0)=l/(1+P!(l)+P!(2)+P!(3)) 

325  P!(1)=P!(1)*P!(0) 

326  P!(2)=P!i2)*P!(0) 

327  P!(3)=P!(3)*P!(0) 

330  IF  NUMRMAN=1  THEN  PD0HN!=P! (2)+P! (3):LQ'=PD0WN'+P! (3)  ELSE  LQ'=P! (3):PD0XN'=P! (3) 

340  L'=P!(1)+2»P!(2)+3*P!(3) 

350  LAMBDA!=(2*P!(0)+2*P!(1)+P!(2))«ALPHA! 

360  MQ!=LQ! /LAMBDA! 

370  W'=L! /LAMBDA' 

380  PSTDBY!=P'(0) 

390  AVL!=P:(01+(P'(l)+P!(l)+P!(2))/3 

395  ' 

499  -=====PRINT  MAIN  DISPLAY  R0UTINE====================== 

500  SCREEN  1,0: COLOR  9,0 

502  IF  FIRSTSCRN  THEN  60SUB  15000: FIRSTSCRN=FALSE: GOTO  510 
505  PUT  (0,0),MAINSCRN 

510  LOCATE  14,27:PRINT  USING  FHT2$; 1/ALPHA!; 

511  LOCATE  14,37:PRINT  USING  FMT2$;l/MUf; 

512  Y=33:G0SUB  10100:  '  Print  statistics  to  screen. 
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520  LOCATE  12,37:PRINT  USING  FHT2$;TIHE; 

521  GOTO  1000 

522  ' 

899  '=====PR06RAH  END  RQUT1NE============================== 

900  SCREEN  O,0:HIDTH  80:C0L0R  6:PRINT  'Proqrai  ended. ..':END 

998  ' 

999  '=====CL0CK  H0DULE====================================== 

1000  FAILTIHEf  (1)=1E+31:  'Calculate  1st  failure  tue  (B  I  0 
1010  DEF  SEG  =iHlAOO:CALL  RNGEN(UI) 

1012  FAILTIHE'(2)=  FNEXPONT! (U! , ALPHA! )+TIME 

1020  CALL  RNGEN(U') 

1022  FAILTIHE'(3)=  FNEXPONT! (U! , ALPHA! 1+TIHE 

1030  REPRTIHE! (1>=1E+31:  '1E+31  =  "infinity" 

1032  REPRTIHE! (2)=1E+31 

1034  REPRTIHE! (3)=1E+31 

1036  STDBYHELQ$="A" 

1099  '**  Determine  next  event  ** 

1100  WHILE  NXTEVNK9950:  'Sets  upper  liiit  to  prgi  run  tin. 
1110  IF  NOT  STATE=0  THEN  1120 

1112  NXTFAIL'=  FNI1IN!(FAILTIHEI(1),FAILTIHE!(2),FAILTIHE,(3)) 

1114  NXTREPR1  =  1E+31 

1116  GOSUB  1021O:BOT0  11B0 

1120  IF  NOT  STATE=3  THEN  1130 

1122  NXTFAIL!=  1E+31 

1124  NXTREPR'=  FNHIN' (REPRTIHE! (1) , REPRTIHE! (2) , REPRTIHE1 (3) ) 

1126  GOSUB  10220:GQT0  1180 

1130  NXTFA1L'=  FNHIN' (FAILTIHE • ( 1 ) , FAILTIME ' (2) ,FAILTI«E' (3) ) 

1132  NXTREPR'=  FNMIN1 (REPRTIHE1 (1) , REPRTIHE! (2) .REPRTIHE! (31) 

1134  GOSUB  10210 

1136  GOSUB  10220 

1180  IF  NXTFAIL!<=MXTREPR'  THEN  1184 

1182  NXTEVNT  =  FNROUNDUP(NXTREPR' 1 :GOT0  1200 

1 184  NXTEVNT  =  FNROUNDUPINXTFAIL!) 

1185  ' 

1199  '**  Run  clock  and  update  tiae  chart  *♦ 

1200  DEF  SEG:P0KE  iH6A,0:  'Clear  key  buffer 
1210  WHILE  T1HE  <  NXTEVNT 

1220  TIHE=TIHE+t 

1230  IF  T1H£<=50  THEN  60SUB  11050: XX=XX+5:F0R  1=1  TO  900: NEXT: GOTO  1240 

1232  DEF  SEG  =  IHIAOO 

1234  CALL  SCRNSHIFT 

1238  XX=276:G0SUB  11050 

1240  LOCATE  12,37:PRINT  USING  FHT2$;TIHE; 

1241  A$=INKEY$  :IF  A$=""  THEN  1260 

1242  IF  A$="s"  OR  A$="S"  THEN  120 

1244  IF  A$="p"  OR  A$="P"  THEN  GOSUB  10010 
1250  DEF  SE6:PQKE  tH6A,0:  'Clear  key  buffer 

1260  WEND 

1261  ' 

1299  '♦*  Update  display  and  determine  next  event  *♦ 

1300  SOUND  120,3:SUMP' (STATE ) =SUMP ! (STATE ) +TIHE-LASTEVNT 
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1302  LASTEVNT  =  NXTEVNT 

1310  IF  NXTFAIL!<=MXTREPR!  THEN  GOSUB  2000  ELSE  GOSUB  3000 

1311  ' 

1399  '«  Update  statistical  estimates  ** 

1400  LAH8DA!=SUMFAIL/TIHE 

1410  IF  NUHRHAN=1  THEN  PDOWN!= (SUMP! (2)+SUMP! (3) ) /TIHE : LQ ! =PDQHN ! +SUMP ! (3) /TIME  ELSE  LQ!=SUMP! (3)/TI 

HE:PDOWN'=LQ! 

1420  L'=(SUMP! (D+SUHP! (2)+SUMP' (2) +SUMP ! (3>  +SUMP » (3>  +SUMP ! (3) > /TIME 

1430  UQ ' =SUMHQ/SUMFAIL: H'=«Q ! +1 /MU ! : PSTDBY '=SUMP ! (0) /TIME 

1440  AVL!=SUMP! (0) /TIME+ (SUMP ! (D+SUHP! ( 1 ) +SUMP ' (2) )/ (T IME+TIME+TIME) 

1450  Y=25:GQSUB  10100:  'Print  estiiates  to  screen 

1451  ' 

1490  A*=  INKEY*: IF  A$="  THEN  1500 

1491  IF  A*="s"  OR  W="S*  THEN  120 

1493  IF  A$="p"  OR  A$='P"  THEN  GOSUB  10010 

1495  DEF  SEG:  POKE  iH6A,0:  'Clear  key  buffer 

1500  WEND  :'No  stop  coaaand  so  go  back  to  line  1100. 

1510  CLS:  LOCATE  12,1:  PRINT  "Prograa  tiae  liait  reached." 

1512  PRINT  "Press  any  key  to  return  to  aenu..." 

1514  GOSUB  10000:  GOTO  120 

1515  ' 

1999  '====FAILURE  H0DULE=======================™*===== 

2000  SUHFA I L=SUHF A IL ♦ 1 : STATE  =  STATE  +  1 

2020  IF  FIRSTFAILURE  THEN  GOSUB  121 00: FIRSTFAILURE=FALSE: GOTO  2130 
2030  BZ=  INSTRCABCFAILINGHELOI) 
2032  X=  FNHELOPOSIT(FAILINGHELOI) 

2040  ON  STATE  GOTO  2100,2200,2300 

2041  " 

2100  DEF  SEG  =&H1A00:CALL  RNGEN(U!) 

2110  REPRTIME* (A2)=  FNEXPONT! (U! ,MU!)+TIHE 

2112  FAILTIME! (AZ)=  1E+31 

2120  GOSUB  12220:  'Ground  flying  helo  i  assign  repairaan. 

2130  AI=  INSTR("A8C",STDBYHEL0$) 

2132  X=  FNHELOPOSIT(STDBYHELO$) 

2140  CALL  RNGENtU!) 

2142  FAILTIME!  IAX)=  FNEXPONT! (U! .ALPHA! 1+TIHE 

2150  GOSUB  12230:  '  Move  stand-by  helo  to  flying  status. 

2160  STDBYHELO$="":RETURN 

2161  ' 

2200  IF  NUMRHAN  =2  THEN  2240 

2210  GOSUB  12210:  'Ground  flying  helo:  no  repainan. 

2220  HAITREPR*  =  HID*(*ABC',AX,1) 

2222  REPRTIME' (AD  =  1E+31 : FAILT IME ' (AI)  =  1E+31 

2224  STARTQUE(AI)  =  TIME 

2230  RETURN 

2240  GOSUB  12220:  'Ground  flying  helo  I  assign  a  repairaan. 

2250  DEF  SEG  =IH1A00:  CALL  RNGENtU!) 

2252  REPRTIME' (AI)  =  FNEXPONT! (U! ,HU! )+TIME 

2254  FAILTIME! (AI)  =  1E+31:  HAITREPR*=" 

2260  RETURN 

2261  ' 
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2300  REPRTIHEMAD  =  I E+3 1 : FAILT I KE ! (AD  =  1E+31 
2310  IF  NUMRHAN=2  THEN  HAITREPR*=  MID* ( "ABC" , Al, 1 ) 
2315  STARTQUE(Al)  =  TIME 
2320  GOSUB  12210:  'Ground  flying  helo:  no  repainan. 

2330  RETURN 

2331  ' 

2999  -=====REPAIR  SUBROUTINE================================ 

3000  AI=  INSTR ( "ABC" ,REPAIREDHELO$) 
3010  X  =  FNHELOPOSIT (REPAIREDHELO$) 

3020  ON  STATE  GOTO  3100,3200,3300 

3021  ' 

3100  STDBYHELO$=  HID$CABC',AI,1) 

3110  FAILTIHE! (AD  =  1E+31 

3120  REPRTIHEMAD  =  1E+31 

3130  60SUB  13010:  'Change  red  helo  to  cold  stand-by  helo. 

3140  STATE  =  STATE-1:  RETURN 

3141  ' 

3200  DEF  SEG  =iHlA00:CALL  RNGEN(U') 

3210  FAILTIHE1 (AD  =  FNEXPQNT' (U! , ALPHA! )+TI«E 

3212  REPRTIHEMAD  =  1E+31 

3220  60SUB  13020:  'Change  red  helo  to  green  and  fly. 

3230  IF  NUHRHAN=2  THEN  STATE=STATE-1:  RETURN 

3240  A2=  INSTRCABC'.HAITREPR*) 

3242  K  =  FNHELOPOSIT (HAITREPR$) 

3244  SUHW9  =  SUrtMQ+TIME-STARTQUE (AD 

3250  GOSUB  13030:  'Change  red  helo  to  red  helo  w/  repairtan. 

3260  CALL  RNGEN(U!  ):REPRTIME' CAD  =  FNEXPONT! (U! ,MU!)+TIME 

3262  WAITREPR$=" 

3270  STATE=STATE-1:  RETURN 

3271  ' 

3300  DEF  SEG  =M1A00:CALL  RNGENfU!) 

3310  FAILTIHE! (AX) =  FNEXPONT' (U! , ALPHA !) +TIHE 

3312  REPRT1ME!  (AD=  1E+31 

3320  GOSUB  13020:  'Change  to  green  helo  and  fly. 

3330  A1X=  INSTRCABC',HAITREPR$) 

3332  K=  FNHELOPOSIT(HAITREPR$):  'Helo  at  head  of  queue. 

3334  SUMWQ  =  SUMWQ+TIME-STARTQUE (AID 

3340  GOSUB  13030:  'Change  red  helo  to  red  helo  w/  repairman. 

3342  CALL  RNGENfU!) 

3344  REPRTIME!  (A1D=  FNEXPONT1  (U!  ,HU')+TINE 

3346  IF  NUMRMAN=2  THEN  WAITREPR$=" : GOTO  3390 

3350  ON  AX  GOTO  3362,3372,3382  :'  Advance  the  queue. 

3360  'REPA1REDHELQS  =  "A" 

3362  IF  HAITREPR$='B"  THEN  KAITREPR*=T  ELSE  «AITREPR$='B' 

3364  GOTO  3390 

3370  'REPAIREDHELO*  =  'B* 

3372  IF  HAITREPR^'A"  THEN  HAITREPR$=T  ELSE  HAITREPR$='A' 

3374  GOTO  3390 

3380  'REPAIREDHELOJ  =  T 

3382  IF  WAITREPR$="A"  THEN  HAITREPR$="B"  ELSE  WA ITREPR$='A' 

3390  STATE=STATE-1:  RETURN 
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3391   ■ 

9998  '*«♦**  SUBROUTINES  *****tmt»t*t»*t**ttt»t*t*t***m* 

9999  ' 

10000  DEF  SEBrPOKE  iH6A,0:  'General  input  routine========== 

10001  A$=  INKEY$:IF  A$="  THEN  10001  ELSE  RETURN 

10002  ' 

10010  DEF  SEG:  POKE  iH6A,0:  'Pause  routine================== 

10011  Al=  INKEYJ:  IF  A$=*c'  OR  M=T   THEN  RETURN  ELSE  GOTO  10011 

10012  ' 

10099  '=====STATISTICS  PRINT  SUBROUTINE=====— ============== 

10100  LOCATE  1B,Y:PRINT  USING  FMTl$; LAMBDA ! ; 

10101  LOCATE  19,Y:PRINT  USING  F«T15;HQ! ; 

10102  LOCATE  20,Y:PRINT  USING  FMT1$;LQ! ; 

10103  LOCATE  21,Y:PR1NT  USING  FMTU;W!; 

10104  LOCATE  22,Y:PRINT  USING  FMT1$;L'; 

10105  LOCATE  23,Y:PRINT  USING  FHT1$;AVL!; 

10106  LOCATE  24,Y:PRINT  USING  FHT1$; PDOWN! ; 

10107  LOCATE  25,Y:PRINT  USING  FMT1$; PSTDBY ! ; 

10108  RETURN 

10109  ' 

10209  '====NEXT  EVENT  SUBR0UT1NES=======================*== 

10210  IF  NXTFAIL!=FAILTIHE!(1)  THEN  FAILINGHELOI^A'  ELSE  IF  NXTFAIL!=FAILTIHE!(21  THEN  FAILINGHELO* 
='B*  ELSE  FAIL1NGHEL0$='C" 

10211  RETURN 

10220  IF  NJTREPR!=REPRTIHE!(n  THEN  REPAIREDHELQ^'A'  ELSE  IF  NXTREPR ! =REPRTI«E ! (2)  THEN  REPAIREDHEL 
0$=,B'  ELSE  REPAIREDHELO*=,C 

10221  RETURN 

10222  ' 

11049  '=====TIC  HARK  SU8R0UT 1  NE========================== 

11050  A2=STATE+1:0N  AI  GOTO  11051,11052,11053,11054 

11051  Y=75:L1NE  (XX, Y)  -  STEP  (5,0) ,2:RETURN 

11052  Y=59:LINE  (XX, Yl  -  STEP  (5,0) ,2:RETURN 

11053  Y=43:LINE  (XX, Y)  -  STEP  (5,0) ,2:RETURN 

11054  Y=27:LINE  (XX, Y)  -  STEP  (5,0) ,2:RETURN 

11055  ' 

12099  '=====FIRST  FAILURE  SU6R0UTINE======================= 

12100  DEF  SEG  =iHlAOO:CALL  RN6EN(U!) :A'=  FNEXPONT' (U1 ,MU')+TIHE 

12102  IF  FAILINGHELO$='B'  THEN  X=54:X1=78:REPRTIHE! (2) =A! : FAILTIHE • (2)=1E+31:G0T0  12106 
12104  IF  FAILINGHELOf="C'  THEN  X=l 07: X 1=1 3 1 : REPRT I  ME ! (3)=A!:FAILTIHE! (3)  =  1E+31 
12106  PUT  (X,105),GRNHELO 

12110  PI  =  3.141593 

12111  CIRCLE  (XI, 161), 12, 2,0, PI 

12112  DRAW  'C2  BU12  U4  M-12,  +  10  M-8,-4  D10  BH+20,-16  H+12,  +  10  R6  fl+4,+5" 

12113  CIRCLE  (Xl,165),14,2,7.5*PI/6,U.5tPI/6 

12114  DRAW  "C2  BD2  L4  NN-2,-15  M-8,-£>" : DRAW  "Dl  1  H-4,+5  D4  NL1  NR1  U3  H+6,-3" 

12115  DRAW  "BR20  H+6,+3  D3  NL1  NR1  U4  H-4,-5  U11":DRAH  'M-8,+6  NH+2,-15  L4" 
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12116  PAINT  STEP  (0,41,2,2 

12117  DRAW  "CO  LA  D6  Rl  U6  R6  D6  Rl  U6  L4" 

12118  CIRCLE  STEP  (0,4) ,2,0:PSET  STEP  (0,0), 0 
12120  GET  (1,143)-  STEP  (48,40) , REDHELO: X 1=X1 +19 

12130  CIRCLE  (X1,165),4,1:PAINT  1X1,1651,3,1 

12131  DRAW  "CI  BM-4,-2  NL3  M+7,+1  M-4,-2' 

12132  LINE  STEP  (0,6)-  STEP  (2,1),3,BF 

12133  LINE  STEP  (1,0)-  STEP  (-4,6),1,BF 

12134  LINE  STEP  (1,0)-  STEP  (2,6),1,BF 

12135  DRAW  "CI  L3  8U9  L2  M-4,-1  Ul  M+4,+1  R2':PAINT  STEP  (-1,1), 1,1 

12140  SET  (X,143)-  STEP  (48,40! , REDHELQ2: RETURN 

12141  ' 

12209  =====FAILURE  DRAWING  SUBROUT I NES====================== 

12210  PUT  a,105),6RNHEL0:PUT  (1,143) , REDHELO:  RETURN: '  No  repainan 

12211  ' 

12220  PUT  IX,105),GRNHEL0:PUT  (X, 143) .REDHELO 

12221  FOR  1=1  TO  300:NEXT:PUT  (X,  143) , REDHELO 

12222  PUT  IX,143) , REDHEL02: RETURN:  '  With  repairman 

12223  • 

12230  PUT  (X,143),GRNHEL0:PUT  (X, 105) .GRNHELO: RETURN: '  Cold  standby  to  flying  status. 

12231  ' 

13009  =====REPAIR  NODULE  SUBROUT INES===================== 

13010  PUT  (X,143),REDHEL02:PUT  (X,1431 ,6RNHEL0: RETURN: "  To  cold  standby 

13011  ' 

13020  PUT  (X,143),REDHEL02:  PUT  (X,143) , GRNHELO 

13021  FOR  1=1  TO  200:NEXT:  PUT  (X, 143) , GRNHELO 

13022  PUT  (X, 105), GRNHELO:  RETURN:  'Repaired  and  goes  flying. 

13023  ' 

13030  FOR  1=1  TO  200:NEXT:PUT  (X, 143) , REDHELO:  PUT  (X,143) ,REDHEL02:  RETURN:  "Assign  a  repainan. 

13031  ' 

14999  =====MAIN  SCREEN  DRAWING  SUBROUT I NE ================== 

15000  LINE  <0,8)-(319,8),2:LlNE  (0,10)-(319,10) ,2 

15002  LINE  (31,23)-(31,80):F0R  Y=27  TO  75  STEP  16:LINE  (27,Y)-(31,Y) :NEXT 

15004  LINE  (3l,80)-(281,80) 

15006  FOR  X=36  TO  281  STEP  5 

1500B  IF  (X=81)  OR  (X=131)  OR  (X=18l)  OR  U=231)  OR  (X=281)  THEN  LINE  (X,80)-(X,85)  ELSE  LINE  (X,80) 

-(X,83) 

15010  NEXT 

15012  LINE  (0,100)-(319,1001,2:L1NE  (0, 102)- (319, 102) ,2 

15014  LINE  (161,102)-(161,199),2 

15016  LOCATE  1,8:PRINT  'MACHINE  -  REPAIRMAN  MODEL' 

15018  LOCATE  4,3:PRINT  "3":L0CATE  6,3:PRINT  '2':L0CATE  8,3:PRINT  'l':LOCATE  10,3:PRINT  '0' 

15020  LOCATE  3,1:PRINT  'N(t) '-.LOCATE  12,32:PRINT  'Tiie:'; 

15022  'Helocopter  drawing  routine 

15024  PI  =  3. 141593: X=25:Y2=165 

15026  GOSUB  1503B 

15028  X=78:Y2=127 
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15030  GOSUB  15038 

15032  X=131 

15034  GOSUB  15038 

15036  GOTO  15062 

15038  Yl=Y2-4 

15040  CIRCLE  (I,Yi),12,l,0,PI 

15042  DRAW  "CI  BU12  U4  L22  R44" 

15044  CIRCLE  (X,Y2) , 14, 1 ,7.5*Pl/6, 1 1.5*PI/A 

15046  DRAW  "CI  BD2  L4  NH-2,-15  H-8,-6" 

15048  DRAW  "Dll  M-4,+5  D4  NL1  NR1  U3  H+6,-3" 

15050  DRAW  "BR20  M+6,+3  D3  NL1  NR1  U4  M-4,-5  Ull" 

15052  DRAW  'M-8,+6  NM+2,-15  L4" 

15054  PAINT  STEP  (0,4), 1,1 

15056  DRAW  "CO  L4  D6  Rl  U6  R6  D6  Rl  U6  L4" 

15058  CIRCLE  STEP  10,4)  ,1,0:PSET  STEP  (0,0) ,0 

15060  RETURN 

15062  'Statistics  Routine 

15064  LOCATE  14,22:PRINT  "HTTF:";: LOCATE  14,32:PRINT  "MTFR:"; 

15066  LINE  (161,113)-(319,113),2 

15068  LOCATE  16,22:PRINT  "Var":LOCATE  16,26:PRINT  "Estiiate^LOCATE  16,35:PRINT  "Luif 

15070  LOCATE  18,22:PRINT  "AR" 

15072  LOCATE  19,22:PRINT  "Hq* 

15074  LOCATE  20,22:PRINT  "Lq" 

15076  LOCATE  21,22:PRINT  "W 

15078  LOCATE  22,22:PRINT  "L* 

15080  LOCATE  23,22:PRINT  "Avl" 

15082  LOCATE  24,4:PRINT  'ALLOCATE  24,10:PRINT  "B"; 

15084  LOCATE  24,17:PRINT  "C";:LOCATE  24,22:PRINT  "Pd"; 

15086  LOCATE  25,22:PRINT  "Ps"; 

15087  GET  (0,0)  -  (319,199) ,HAINSCRN 

15088  GET  (1,143)-(49,183),GRNHEL0 
150B9  RETURN 

15090  ' 

19999  -=====TITLE  SCREEN================================= 

20000  SCREEN  0,1:  COLOR  14,1,1:  WIDTH  40:  CLS 
20010  LOCATE  2,12:  PRINT  'GRAPHIC  SIMULATION* 
20012  LOCATE  4,18:  PRINT  "of  the' 

20014  LOCATE  6,9  :  PRINT  -MACHINE  -  REPAIRMAN  MODEL" 
20016  LOCATE  9,14:  PRINT  'by  R.E.  Nelsen" 

20020  LOCATE  14,5:  PRINT  "Submitted  in  partial  fulfillment* 

20021  LOCATE  15,2:  PRINT  "of  the  requireients  for  the  degree  of 

20022  LOCATE  17,1:  PRINT  "Master  of  Science  in  Operations  Research" 

20023  LOCATE  19,4:  PRINT  "froi  the  Naval  Postgraduate  School" 

20024  LOCATE  20,11:  PRINT  "Monterey,  California" 

20025  LOCATE  22,2:  PRINT  "Advisors:  J.D.  Esary,  A.F.  Andrus" 

20030  LOCATE  25,2:  PRINT  "Press  any  key  to  continue... ";:GOSUB  10000 

20040  CLS:LOCATE  1,14:  PRINT  "Acknowledgment" 

20041  LOCATE  5,6:PRINT  "The  author  would  like  to  thank" 

20042  LOCATE  7,2:PRINT  "Associate  Professor  Bruno  Shubert  for" 

20043  LOCATE  9,2:PRINT  "the  donation  of  his  asseibly-language" 
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20044  LOCATE  11,2:PRINT  'random-number  generator  which  formed" 

20045  LOCATE  13,2:PR1NT  "the  basis  for  RHGEN.SRT." 

20046  LOCATE  25,2:PRINT  'Press  any  key  to  continue. .."; :GQSUB  10000:RETURN 

20047  ' 

20099  '  =====  I NSTRUCT 1 0NS================================= 

20100  SCREEN  0,1: COLOR  14, 1 , IsHIDTH  40:CLS 

20110  LOCATE  1,4:PRINT  'Instructions  lor  the  Prograi  Menu'; 

20111  LOCATE  3,2:PRINT  "Press  <s>  to  set  the  random-number  seed"; :PRINT  "and  enter  a  positive  intege 
r  when  asked. The  prograi  returns  to  the  ■enu." 

20112  PRINT:PRINT  '  Press  <c)  to  change  model  parameters.  The  prograi  will  ask  for  the  machine's  i 
ean  tiie  to  failure  (HTTF),  the  repair-"; 

20113  PRINT  "ten's  mean  tile  for  repair  (HTFR),  and  the  number  of  repairmen  to  employ.  The" 

20114  PRINT  "HTTF  and  the  HTFR  should  be  positive   integers  that  are  greater  than  10.  The  number 

of  repairmen  can  be  1  or  2.  The  simulation  starts  after  the  number  of   repairmen  has  been  entered. 

■ 

20115  LOCATE  18,2:  PRINT  'Press  <d>  to  use  default  parameters. ";:PRINT  'The  simulation  will  immediat 
ely  start   with  the  values:  HTTF  =  15,  HTFR  =  10,  number  of  repairmen  =  1." 

20116  LOCATE  23,2:  PRINT  "To  end  the  program,  press  <e>.' 

20117  LOCATE  25,2:  PRINT  'Press  any  key  to  continue. .. ";:G0SUB  10000 

20120  CLS:L0CATE  1,13:  PRINT  "Keyboard  Commands'; 

20121  LOCATE  4,2:  PRINT  "During  the  simulation,  you  may  use  the";:PRINT  "following  keyboard  commands 

.  ■ 

20122  LOCATE  7,5:  PRINT  '<p>  -  pause  the  simulation.' 

20123  LOCATE  9,5:  PRINT  "<c>  -  continue  the  simulation." 

20124  LOCATE  11,5:  PRINT  "<s>  -  stop  the  simulation  and  return" 

20125  LOCATE  13,11:  PRINT  'to  the  program  menu." 

20126  LOCATE  25,2:PRINT  "Press  any  key  to  continue... " ; : G0SUB  10000 

20130  CLS:L0CATE  1,13:  PRINT  'Graphic  Display"; 

20131  LOCATE  3,2:PRINT  'The  variable  names  in  the  display  are";:PRlNT  'defined  as:"; 

20132  LOCATE  6,1:  PRINT  'TIHE:  current  time  on  simulation  clock."; 

20133  LOCATE  8,1:  PRINT  *N(t):  no.  of  machines  down  at  time  t."; 

20134  LOCATE  10,3:  PRINT  "AR:  average  arrival  rate  of  machines"; 

20135  LOCATE  11,7:  PRINT  "into  the  repair  queue."; 

20136  LOCATE  13,3:  PRINT  'Hq:  average  waiting  time  of  machines'; 

20137  LOCATE  14,7:  PRINT  'in  the  repair  queue."; 

20138  LOCATE  16,3:  PRINT  "Lq:  average  length  of  repair  queue."; 

20139  LOCATE  18,4:  PRINT  "W:  average  down  time  of  a  machine.'; 

20140  LOCATE  20,4:  PRINT  'L:  average  no.  of  down  machines. "; 

20141  LOCATE  22,2:  PRINT  'Avl:  machine  availability.'; 

20142  LOCATE  25,2:  PRINT  'Press  any  key  to  continue... ";:60SUB  10000 

20143  CLS:L0CATE  1,3:  PRINT  "Pd:  probability  that  a  down  machine"; 

20144  LOCATE  2,7:  PRINT  'must  wait  for  repair.'; 

20145  LOCATE  4,3:  PRINT  'Ps:  probilility  that  an  up  machine"; 

20146  LOCATE  5,7:  PRINT  "must  wait  in  cold  standby."; 
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20147  LOCATE  8,1:  PRINT  '  The  nuabers  listed  under  the  coluan  heading  'Liait'  are  the  halting  va 
lues  as  tiie  goes  to  infinity  of  the  above  defined  quantities.  These  values  ire  computed  with  t 
he  balance  equations  for  the  aodel." 

20148  LOCATE  15, 1: PRINT  '  The  values  listed  under  the  CDluan  heading  'Estiaate'  are  the  estimate 
s  of  the  liaiting  values.  These  quantities  are  calculated  with  data  collected  frot' 

20149  PRINT  "the  siaulation.  The  estnates  are  recoa-puted  upon  every  occurance  of  an  event." 

20150  LOCATE  25,2:PRINT  "Press  any  key  to  return  to  ■enu...,;:GQSUB  10000:GOTO  120 

20151  ' 

29999  '=====BASICA  AND  COLOR/GRAPHICS  ADAPTER  CHECK========= 

30000  DEF  SEG  =0: IF  (PEEKI&H410)  AND  IH30)  <>  iH30  THEN  DEF  SEG:G0T0  30007 

30003  LOCATE  3,1:PRIMT  "Sorry..." 

30004  PRINT  "You  do  not  have  the  color/graphics  aonitor  adapter!" 

30005  PRINT  "This  siaulation  uses  graphics  and  requires  that  adapter." 

30006  DEF  SEG  :  END 

30007  ON  ERROR  GOTO  3000B:  PLAY  "pl6":GOT0  30011 

30008  WIDTH  80:CLS:L0CATE  3,1 

30009  PRINT  "This  siaulation  uses  advanced  BASIC." 

30010  PRINT  'Reload  this  prograa  after  using  the  coaaand  'BASICA'. ":END 

30011  ON  ERROR  GOTO  0:  RETURN 

30012  ' 

30019  '===~FILE-L0ADING  ERROR  TRAPPING  R0UTINE============= 

30020  IF  ERR  <>  53  THEN  ON  ERROR  GOTO  0 

30022  CLS:LOCATE  3,1:PRINT  "On  which  drive  can  this  prograa  be  found?" :G0SUB  10000 
30024  DRVE$=A$+':":  RESUME  82 
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APPENDIX  D 


RNGEN.SRT  SOURCE-CODE  LISTING 


0000 


0000 

0000 

0000 

55 

0001 

8B  EC 

0003 

IE 

0004 

8C  C8 

0006 

8E  DB 

0008 

Al  0161  R 

000B 

8B  16  0163  R 

000F 

8B  F0 

0011 

SB  FA 

0013 

83  FA  01 

0016 

77  6A 

0018 

72  05 

001A 

3D  F31D 

001D 

73  63 

001F 

Dl  EO 

0021 

Dl  D2 

0023 

Dl  EO 

0025 

Dl  D2 

0027 

Dl  EO 

0029 

Dl  D2 

imnmiiniuiniiiiiiiiiiuiiniiiiniiiiinn 
RNGEN.SRT  -  BASIC  USR  Routine 
by  B.O.  Shubert 
Modified  by  R.E.  Nelsen 
2  flay  1984 

This  subroutine  generates  and  returns  a  unifors 
randot-nuaber  in  BASIC'S  single-precision  foraat. 
The  seed  X  is  not  uported  nor  returned  to  BASIC, 
but  it  is  autoaatically  updated.  The  algontha 
used  is: 

X  =  AX  nod  M 

with  A  =  7A5  and  H  =  2A31  -  1.  RNGEN  is  called 
froa  BASIC  using  the  statements; 

RNGEN  =  0 

DEF  SEG  =  IH1A00 

CALL  RNGEN  (U1) 
U!  is  the  returned  randoa  nuaber  and  aust  be 
defined  prior  to  the  first  call  to  RNGEN. 
ininiiiiiijiiiiiinmimiiiiiiiiiuiiimuni 


RNGEN 


MULT32: 


CSEG 

SEGMENT 

ASSUME  CS : CSEG , DS : CSES , ES : NOTH I NG 

0R6 

0 

PROC 

FAR 

PUSH 

BP 

;save  for  BASIC 

HOV 

BP,SP 

jpoint  to  args  on  stack 

PUSH 

DS 

;save  for  BASIC 

NOV 

AX,CS 

;aake  data  accessable 

HOV 

DS,AX 

NOV 

AX,L0SEED 

;get  old  X 

HOV 

DX.HISEED 

NOV 

SI ,  AX 

nov 

DI.DX 

CHP 

DX,  1 

;check  if  AIX  <  M 

JA 

HULT48 

;48  bit  aultiply  if  not 

JB 

HULT32 

;else  32  bit  aultiply 

CUP 

AX.0F31DH 

JNB 

HULT48 

SAL 

AX,  1 

;aultiply  X  by  7A5 

RCL 

DX,  1 

SAL 

AX,1 

RCL 

DX,  1 

SAL 

AX,  1 

RCl 

DX,  1 

78 


002B 

2B  C6 

SUB 

AX, SI 

002D 

IB  D7 

SBB 

DX.DI 

002F 

8B  F0 

MOV 

SI, A? 

0031 

BB  FA 

MOV 

DI,DX 

0033 

Dl  EO 

SAL 

AX,  1 

0035 

Dl  D2 

RCL 

DX,  1 

0037 

Dl  EO 

SAL 

AX,  1 

0039 

Dl  D2 

RCL 

DX,  1 

003B 

Dl  EO 

SAL 

AX,  1 

003D 

Dl  D2 

RCL 

DX,  1 

003F 

2B  C6 

SUB 

AX, SI 

• 

0041 

IB  D7 

SBB 

DX,DI 

0043 

8B  FO 

HOV 

SI, AX 

0045 

8B  FA 

HOV 

DI,DX 

0047 

Dl  EO 

SAL 

AX,  1 

0049 

Dl  D2 

RCL 

DX,  1 

004B 

Dl  EO 

SAL 

AX,  1 

004D 

Dl  D2 

RCL 

DX,  1 

004F 

Dl  EO 

SAL 

AX,  1 

0051 

Dl  D2 

RCL 

DX,  1 

0053 

2B  C6 

SUB 

AX, SI 

0055 

IB  D7 

SBB 

DX,DI 

0057 

8B  FO 

NOV 

SI, AX 

0059 

BB  FA 

MOV 

DI,DX 

005B 

Dl  EO 

SAL 

AX,  1 

005D 

Dl  D2 

RCL 

DX,  1 

005F 

Dl  EO 

SAL 

AX,  1 

0061 

Dl  D2 

RCL 

DX,  1 

0063 

Dl  EO 

SAL 

AX,  1 

0065 

Dl  D2 

RCL 

DX,  1 

0067 

2B  C6 

SUB 

AX, SI 

0069 

IB  D7 

SBB 

DX.DI 

006B 

BB  FO 

HOV 

SI,  AX 

006D 

BB  FA 

HOV 

DI,DX 

006F 

Dl  EO 

SAL 

AX,  1 

0071 

Dl  D2 

RCL 

DX,  1 

0073 

Dl  EO 

SAL 

AX,  1 

0075 

Dl  D2 

RCL 

DX,  1 

0077 

Dl  EO 

SAL 

AX,  1 

0079 

Dl  D2 

RCL 

DX,  1 

007B 

2B  C6 

SUB 

AX, SI 

007D 

IB  D7 

SBB 

DX,DI 

007F 

E9  012F  R 

JHP 

UNIFOi 

;done  since  AtX^H 

0082 

33  DB 

! 

HULT48:   XOR 

BX,BX 

jwltiply  X  by  7A5 

0084 

8B  CB 

HOV 

CX,BX 

0086 

Dl  EO 

SAL 

AX,  2 

;by  multiplying  X  by 

0088 

Dl  D2 

RCL 

DX,  1 

008A 

Dl  D3 

RCL 

BX.l 
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008C 

Dl  EO 

OOBE 

Dl  D2 

0090 

Dl  D3 

0092 

Dl  EO 

0094 

Dl  D2 

0096 

Dl  D3 

0098 

2B  C6 

009A 

IB  D7 

009C 

IB  D9 

009E 

8B  FO 

OOAO 

8B  FA 

00A2 

SB  CB 

00A4 

Dl  EO 

0OA4 

Dl  D2 

00A8 

Dl  D3 

OOAA 

Dl  EO 

OOAC 

Dl  D2 

OOAE 

Dl  D3 

OOBO 

Dl  EO 

00B2 

Dl  D2 

00B4 

Dl  D3 

00B6 

2B  C6 

00B8 

IB  D7 

OOBA 

IB  D9 

OOBC 

8B  FO 

OOBE 

8B  FA 

OOfO 

8B  CB 

00C2 

Dl  EO 

00C4 

Dl  D2 

00C6 

Dl  D3 

00C8 

Dl  EO 

OOCA 

Dl  D2 

OOCC 

Dl  D3 

OOCE 

Dl  EO 

OODO 

Dl  D2 

00D2 

Dl  D3 

00D4 

2B  C6 

00D6 

IB  D7 

00D8 

IB  D9 

OODA 

SB  FO 

OODC 

SB  FA 

OODE 

8B  CB 

OOEO 

Dl  EO 

00E2 

Dl  D2 

00E4 

Dl  03 

00E6 

01  EO 

OOEB 

01  D2 

OOEA 

Dl  03 

OOEC 

01  EO 

SAL 

AX,  1 

RCL 

DX,  1 

RCL 

BX,1 

SAL 

AX,  1 

RCL 

DX,  1 

RCL 

BX,1 

SUB 

AX, SI 

; subtract  original  X 

SBB 

DX,D1 

SBB 

BX,CX 

MOV 

SI, AX 

jsavinq  ne«  X 

nov 

DI.DX 

• 

HOV 

CX,BX 

SAL 

AX ,  1 

;and  repeat  this  5  tues 

RCL 

DX,  1 

RCL 

EX ,  1 

SAL 

AX,  1 

RCL 

DX,  1 

RCL 

BX,1 

SAL 

AX,  1 

RCL 

DX,  i 

RCL 

BX,1 

SUB 

AX, SI 

SBB 

DX,DI 

SBB 

BX,CX 

HOV 

SI, AX 

NOV 

Dl,  DX 

HOV 

CX,BX 

SAL 

AX,  1 

RCL 

DX,  1 

RCL 

BX,1 

SAL 

AX,  1 

RCL 

DX,  1 

RCL 

BX,1 

SAL 

AX,  1 

RCL 

DX,  1 

RCL 

BX,  1 

SUB 

AX, SI 

SBB 

DX ,  Dl 

SBB 

BX,CX 

nov 

SI, AX 

nov 

Dl ,  DX 

nov 

CX.BX 

SAL 

AX.l 

RCL 

DX,  1 

RCL 

BX,  1 

SAL 

AX,  1 

RCL 

DX,  1 

RCL 

BX,  1 

SAL 

AX,  1 
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OGEE 

Dl  D2 

OOFO 

Dl  D3 

00F2 

2B  C6 

00F4 

IB  D7 

00F6 

IB  D9 

00F8 

8B  FO 

OOFA 

8B  FA 

OOFC 

8B  CB 

OOFE 

Dl  EO 

OlOO 

Dl  D2 

0102 

Dl  D3 

0104 

Dl  EO 

0106 

Dl  D2 

0108 

Dl  D3 

010A 

Dl  EO 

010C 

Dl  D2 

010E 

Dl  D3 

0110 

2B  C6 

0112 

IB  D7 

0114 

IB  D9 

0116 

Dl  E2 

0118 

Dl  D3 

01 1 A 

Dl  EA 

one 

03  C3 

011E 

83  D2  00 

0121 

81  FA  7FFF 

0125 

77  03 

0127 

EB  06  90 

012A 

40 

012B 

81  DA  8000 

012F 

A3  0161  R 

0132 

89  16  0163  R 

0136 

Dl  EO 

0138 

Dl  D2 

013A 

32  C9 

013C 

Dl  EO 

013E 

Dl  D2 

0140 

72  04 

0142 

FE  CI 

0144 

EB  F6 

0146 

Dl  EA 

0148 

Dl  D8 

014A 

B3  80 

014C 

2A  D9 

SUBTRM: 


UNIF01: 


FNDEXP: 


EXPONT: 


RCL 

DX,  1 

RCL 

BX,  1 

SUB 

AX, SI 

SBB 

DX,  Dl 

SBB 

BX,CX 

NOV 

SI, AX 

nov 

Dl,  DX 

HOV 

CX,BX 

SAL 

AX,  1 

RCL 

DX,  1 

RCL 

BX,  1 

SAL 

AX,  1 

RCL 

DX,  1 

RCL 

BX,  1 

SAL 

AX,  1 

RCL 

DX,  1 

RCL 

BX,  1 

SUB 

AX, SI 

SBB 

DX,  Dl 

SBB 

BX,CX 

;BX,DX,AX,  now  contain  AIX 

SAL 

DX,  1 

RCL 

BX,  1 

;get  K=AJX/2A31  to  BX 

SHR 

DX,  1 

;and  X*=ASX (aod  2A31)  to  DX, AX 

ADD 

AX,BX 

calculate  X"  =  X'+K 

ADC 

DX,0 

CttP 

DX,7FFFH 

;see  if  X"  (  B 

JA 

SUBTRM 

JHP 

UNIF01 

;done  if  so  since  X"  can't  =  11 

INC 

AX 

;else  subtract  11  by  adding  1 

SBB 

DX.8000H 

;and  subtracting  2A31 

nov 

LOSEED.AX 

;store  new  X 

nov 

HISEED,DX 

SAL 

AX,  1 

;Divide  X  by  2'31  to  lake 

RCL 

DX,  1 

;a  uniform  (0,1). 

XOR 

CL,CL 

;Set  a  counter  to  zero. 

SAL 

AX,  1 

;Shift  left  and  count  the 

RCL 

DX,  1 

;number  of  leading  zeros. 

JC 

EXPONT 

;Leading  1  found. 

INC 

CL 

;Not  found,  add  1  to  counter 

JHP 

FNDEXP 

;and  look  at  next  bit. 

SHR 

DX,  1 

;Shift  back  to  correct  value 

RCP 

AX,  1 

;and  suppress  leading  1. 

MOV 

BL,80H 

; Compute  normalized  exponent. 

SUB 

BL.CL 

BL,  DX,  and  AX  now  contain  a  single  precision,  uniform 
(0,1)  randoi  number  in  BASIC'S  floating  point  format. 


81 


DUE 

IF 

014F 

BB  7E  06 

0152 

88  25 

0154 

88  55  01 

0157 

88  75  02 

015A 

88  5D  03 

015D 

5D 

015E 

CA  0002 

0161 

4130 

0163 

OOAB 

0165 

0165 

POP 

DS 

jrestore  segment  register 

MOV 

di.ebp: 

1+6 

;get  addr  of  variable  U 

MOV 

[DI],AF 

1 

jpass  LSB  of  aantissa 

nov 

[DIM, 

DL 

jpass  MSB  of  aantissa 

nov 

EDI3+2, 

DH 

;pass  HSB  of  aantissa 

nov 

[DH+3, 

BL 

;pass  exponent 

POP 

BP 

;restore  BP  for  BASIC 

RET 

2 

;FAR  return  to  BASIC 

! 

LOSEED 

DM 

4130H 

;storage  for  K 

HISEED 

DW 

OABH 

RNBEN 

ENDP 
CSEB 
END 

ENDS 
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APPENDIX  E 
SCRNSHFT.SRT  SOURCE-CODE  LISTING 


0000 


0000 

0000 

FD 

0001 

0000 

0003 

0000 

0005 

0086 

0007 

0007 

55 

0008 

FB 

0009 

B9  0024 

OOOC 

BA  001B 

000F 

B4  OD 

0011 

CD  10 

0013 

83  E9  05 

0016 

B4  OC 

0018 

CD  10 

001A 

83  CI  05 

001D 

83  C2  10 

0020 

B4  OD 

0022 

CD  10 

0024 

83  E9  05 

0027 

B4  OC 

0029 

CD  10 

002B 

83  CI  05 

*tttt»tttt{ti»*»fiftt»it<»t*»tf»ttf»*tttfttt* 

SCRNSHFT.SRT,  Ver  2  --  BASIC  USR  Routine 
by  R.E.  Nelsen  --  27  July  1984 

This  subroutine  shifts  5  coluins  to  the  left 
the  4  rows  of  pixels  corresponding  to  the 
■odel's  state  vs  tiie  graph.  It  is  called 
froi  BASIC  using  the  coMands: 

SCRNSHFT  =  0 
DEF  SEG  =  &H1A00 
CALL  SCRNSHFT 

This  subroutines  calls  the  ROM  BIOS  Type  10H 
interrupt  (Video  I/O)  using  routines  12  i  13. 
miitHiiHiiiiiiummitmnmimtnn 

CSEG  SEGMENT 

ASSUME  CS:CSEG,  DS:NOTHING 

Establish  a  file  header  for  use  by  the 
BASIC  BLOAD  coaaand. 

HEADER: 


DB 

OFDH 

;Code  for  BLOAD  file 

DM 

0 

;Seg  addr  location 

DH 

0 

■.Offset  location 

DH 

RTN.LEN 

;Routine  length 

5 

SCRN.SHFT  PROC 

FAR 

PUSH 

BF 

;Save  for  BASIC 

ST  I 

;Enable  interrupts 

MOV 

CX  ,36 

;Scrn  coluin  36 

A6AIN:    MOV 

DX,27 

;Scrn  row  27 

MOV 

AH, 13 

;Read  pix  at  ton  DX 

INT 

10H 

;«  col  CX.  Pix  in  AL 

SUB 

CX,5 

;Shift  left  5  col's 

MOV 

AH, 12 

;Set  =  12  to  write 

INT 

10H 

jpixei  to  new  pos'n 

ADD 

CX,5 

{Shift  right  5  col's 

ADD 

DX ,  16 

;Row  43 

MOV 

AH,  13 

;Set  =  13  to  read 

INT 

10H 

;Read  pixel 

SUB 

CX,5 

;Shift  left  5 

MOV 

AH, 12 

;Set  =  12 

INT 

10H 

;Write  pixel 

ADD 

CX,5 

; Shi  ft  right  5  col "s 
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;Now  repeat  procedure  for  retaining  two  rows. 


002E 

83  C2  10 

0031 

B4  0D 

0033 

CD  10 

0035 

83  E9  05 

0038 

B4  OC 

003A 

CD  10 

003C 

83  CI  05 

003F 

83  C2  10 

0042 

B4  OD 

0044 

CD  10 

0046 

83  E9  05 

0049 

B4  OC 

004B 

CD  10 

004D 

83  CI  06 

0050 

81  F9  01 1A 

0054 

72  66 

ADD 

DI.16 

;Row  59 

HOV 

AH, 13 

INT 

10H 

SUB 

n,5 

MOV 

AH, 12 

INT 

10H 

ADD 

CI, 5 

ADD 

DX.16 

;Row  75 

MOV 

AH, 13 

INT 

10H 

SUB 

CI, 5 

HOV 

AH, 12 

INT 

10H 

ADD 

CI, 6 

;Shift  to  new  coluan 

CMP 

CI, 282 

;Is  last  col  reached? 

JB 

A6AIN 

; If  no,  do  next  col. 

; Fini shed  with  rows,  so  now  iust  blank  oat  last 
;5  colutns  of  each  row. 


0056 

41 

INC 

CI 

;Col  283 

0057 

B4  OD 

HOV 

AH,13 

; Set  to  read  dot 

0059 

CD  10 

INT 

10H 

;Read  blank  dot 

005B 

8A  D8 

HOV 

BL,AL 

;Save  color  code 

005D 

B9  0115 

HOV 

CI, 277 

;Col  277 

0060 

BA  001B 

NEH_CQL:   HOV 

DI,27 

;Row  27 

0063 

B4  OC 

-  HOV 

AH, 12 

;Set  to  write  dot 

0065 

8A  C3 

HOV 

AL,BL 

; Set  backgrnd  color 

0067 

CD  10 

INT 

10H 

; Write  blank  dot 

0069 

B3  C2  10 

ADD 

DI,16 

;Do  again  for 

006C 

B4  OC 

HOV 

AH, 12 

;row  43, 

006E 

8A  C3 

MOV 

AL.BL 

0070 

CD  10 

INT 

10H 

0072 

83  C2  10 

ADD 

DI,16 

;and  row  59, 

0075 

B4  OC 

HOV 

AH, 12 

0077 

8A  C3 

HOV 

AL,BL 

0079 

CD  10 

INT 

10H 

007B 

83  C2  10 

ADD 

DI.16 

;and  row  75. 

007E 

B4  OC 

HOV 

AH, 12 

0080 

8A  C3 

HOV 

AL,BL 

0082 

CD  10 

INT 

10H 

0084 

41 

INC 

CI 

;New  col  on  right 

0085 

81  F9  01 1A 

CHP 

CI, 282 

; Last  coluin7 

0089 

72  D5 

JB 

NEH.COL 

;No,  so  do  next  col. 

» 

;Done,  so  recover  BP  and  return  to  BASIC. 
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008B  5D 
008C  CE 
008D 
=  0086 
008D  1A 
008E  . 


POP 

BP 

RET 

SCRNJHFT 

ENDP 

RTN_LEN 

ECU 

$  -  SCRN_SHFT 

DB 

01AH    ;Need  for  BLDAD 

CSEG 

ENDS 

END 

HEADER 
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